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ALGORITHMIC DESIGN OF PEPTIDES FOR BINDING AND/OR 
MODULATION OF THE FUNCTIONS OF RECEPTORS AND/OR OTHER 

PROTEINS 

5 This application is a continuation-in-part of copending application United States 

Serial No. 09/490,701, filed January 24, 2000, which is incorporated herein in its entirety. 

FIELD OF THE INVENTION 
The invention relates generally to peptide molecules and to methods of designing 
1 0 peptides or peptide-like molecules. More particularly, the invention relates to novel, 
short peptides or peptide-like molecules which have a high probability of binding to 
and/or otherwise modulating the function of polypeptides or proteins, and to methods for 
designing such peptides or peptide-like molecules. 

15 BACKGROUND OF THE INVENTION 

All protein sequences, whether peptides, polypeptides, or proteins, are composed 
of a linear sequence of amino acids joined by peptide bonds. There are twenty naturally 
occurring amino acids, each bearing a chemically unique side chain. Determinants of 
polypeptide interactions, such as those between peptide segments in protein folding or 

20 between protein monomers, are encoded in the one-dimensional sequence of these twenty 
amino acid side chains. For purposes of this application, "peptides" are generally 
considered to be amino acid polymers of not more than 25 amino acids in length; 
"polypeptides" are generally considered to be polymers of between 25 and 50 amino 
acids; and "proteins" are generally considered to be polymers containing more than 50 
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amino acids. One of ordinary skill in the art would appreciate that some overlap among 
these ranges is expected, and minor deviations from these ranges does not in any way 
diminish the scope of the invention. The ''naturally occurring amino acids" are those that 
are encoded for in the genetic code, and which are generally considered to be those found 
5 in all living species to date. 

Net differences in the cumulative energetic contributions of several types of weak 
bonding mechanisms, totaling as little as AG = 5-10 kcal/mol, determine selection and 
stabilization among conformations observed in protein folding, protein-protein 
interactions and the initial phases of substrate-enzyme and ligand-membrane receptor 

1 0 association. In particular, the minimization of AG through the formation of four general 
types of weak bonding mechanisms between amino acid side chains, in the range of AG = 
2-7 kcal/mol, determines the arrangement of protein sequences in three-dimensional 
space, as well as the relative orientations of protein chain aggregates, in aqueous 
environments and at physiological temperatures. The thermal instability of the 

1 5 conformations supported by these low AG, reversible, weak-bonding mechanisms permits 
uncatalyzed, fast searches of configuration space for functionally optimal cooperative 
arrangements within and between polypeptide and protein monomers. The variety of 
weak bond capacities afforded by amino acid side chains determines the range of the 
amino acid sequences' physicochemical property transformations listed in this invention. 

20 The weak bonds ordering polypeptides and proteins in three-dimensional space 

include hydrogen bonds, such as the main chain amino acid carbonyl and imino groups, 
which configure the right-turning a-helices and the parallel and antiparallel p-sheets. 
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They also include the hydrogen and ionic bonds between amino acid side chains, such as 
the hydroxyl groups of serine and threonine, the acidic carboxyl groups of aspartate and 
glutamate, and the basic groups of lysine and arginine. In addition to being distinct with 
respect to the chemical group, these weak hydrogen and ionic bonding influences are also 
5 directionally specific, with bonding angles greater than 30° reducing their influence to 
negligible levels. 

A third but nondirectional type of weak bonding interaction, induced by 
fluctuating charges within a distance of 1-3 A, is called van der Waal forces. These 
;i interactions vary with the size and the extent of mutual geometric fit, but are in the range 

H 10 of 1 -2 kcal/mol. These forces are barely greater than those due to the heat of molecular 
:J motion at room temperature (AG 2 0.6-1 .0 kcal/mol). However, in the specific cases of 

some antibody/antigen interactions and MHC protein/peptide interactions, which involve 
ri water-releasing tight fits between corresponding moieties in suitably shaped binding 

r y pockets, the AGs associated with van der Waals interactions have been estimated to be as 

!: :: I 

P 15 high as 30 kcal/mol. 

A fourth weak bonding mechanism, and the most energetically dominant force on 
three-dimensional polypeptide structure and protein-protein interactions, is termed the 
hydrophobic effect. The hydrophobic effect arises from the much stronger attraction that 
water molecules have for each other than for hydrocarbon groups or molecules. Each 
20 tetrahedrally-coordinated water molecule participates in strong, hydrogen-bonded, 

dipole/dipole interactions with other water molecules that are manifested in the properties 
of water such as its high surface tension, high latent heat and high boiling point. These 
physicochemical features of water molecules afford a large variety of possible atomic 
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arrangements of water (as seen in the large number of different ice types) that in turn 
permit maximizing the entropy and minimizing the free energy of the aqueous solution. 
Spatially distributed (nondirectional) deformations in these hydrogen-bonded 
arrangements of water result from the intrusion of nonpolar, hydrophobic solutes. The 
5 introduction of such molecules into an aqueous solution results in the formation of 
volume-expanding hydration shells composed of hydrogen-bonded cages of multiple 
molecular layers of water ("clathrate structures") around these molecules, in a process 
called "hydrophobic hydration". In aqueous solutions, such deformations in water 
structure are energetically unfavored. For example, the side chains of alanine, valine, 

10 leucine and isoleucine are without effective dipole moments, and therefore cannot 

participate in charge-mediated or hydrogen-bonding interactions with water. As a result, 
these side chains intrude into the aqueous solvent and disrupt the ordered structure of the 
aqueous solvent, resulting in an increase in the overall AG. Amino acids with polar but 
uncharged side chains, such as serine and threonine, may hydrogen bond with a molecule 

1 5 of water, but otherwise undergo the same kind of hydrophobic hydration as the non-polar 
side chains. In the case of amino acids with side chains containing charged groups, such 
as glutamate or lysine, the electrostatic fields associated with these side groups are 
screened by water molecules, such that in an aqueous solution hydrophobic hydration is 
still a prominent characteristic of these amino acids as well. The nonlocal, cooperative 

20 interactions of the hydrogen bonds of the aqueous solvent surrounding these amino acids 
drive the in-line, surface-minimizing attraction between the coherent hydrophobic-phase 
patches of amino acid side chains, thereby maximizing the entropy, and minimizing the 
free energy, of the overall aqueous solution. 
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The importance of the sequential arrangements of amino acid side chain 
hydrophobicities in the determination of peptide and protein secondary structures has 
been established knowledge in protein biology for many decades. The ready availability 
of water for compensatory weak bonding implies that relatively small changes in AG 
5 occur when internal peptide backbone-related, carbonyl-imino hydrogen bonding or side 
chain polar groups are not satisfied. This contrasts with the much greater alteration in AG 
associated with loss of internal hydrophobic bonding, which cannot be compensated by 
the hydrophobically disrupted, aqueous environment. Minimization of hydrophobic free 
energy, AGh P , by water interface-reducing aggregation of nonpolar, hydrophobic amino 

1 0 acid side chain groups adds to the AG of binding that can, collectively, be orders of 

magnitude larger than that predicted by van der Waals theory. Mutually attractive forces 
mediated by hydrophobic surface minimization have been measured by atomic force 
spectroscopy to extend to as great a distance as 60 A, the length scale of synaptic gaps. 
These attractive forces decay less than exponentially with distance. The contribution to 

1 5 the energy of stabilization of the three-dimensional, tertiary structure of protein by AGh P 
minimization due to aggregation of hydrophobic amino side chains has been estimated to 
be in the range of 70%. 

Complete substitution of hydrophobically equivalent amino acids in peptides 
maintains and sometimes increments their peptide-receptor mediated physiological 

20 potency. Additionally, proteins which are dominated by helical secondary structures of 
specific turn lengths can be designed using sequences of amino acids of high and low 
hydrophobicities, independent of the specific amino acids chosen within each 
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hydrophobicity class. In contrast, regions of amino acids characterized by interactions 
dominated by hydrogen bonds, ionic bonds, and van der Waals interactions are often 
exquisitely sensitive to any substitution, even those deemed to be conservative 
replacements. This difference between the effects on AG of hydrophobic interactions 
5 versus those of hydrogen bonding, ionic binding or van der Waals interactions, along 
with more stringent geometric requirements of the latter compared with hydrophobic 
weak bonds, make sequential patterns of AGhpin polypeptide sequences of primary 
importance in determining peptide-peptide or peptide-protein interactions. 

Previously, the role of the hydrophobic interactions of amino acids in peptide 

1 0 ligands with amino acids in their associated membrane proteins have been considered in 
structure-function analyses in two ways. First, the local roles of amino acids have been 
evaluated. In these studies, ligand-receptor binding is changed by point mutations in 
specifically positioned amino acids, producing alterations in the hydrophobic 
characteristics of "binding pockets" involving neighboring but nonsequential 

1 5 juxtapositions of residues brought together in the protein's cooperative tertiary structure. 
Second, the global effects of amino acids have been examined. These effects are often 
studied using chimeric exchanges, with respect to the number, lengths, and locations of 
transmembrane segments of receptors, transporters, and/or channels, and exploit the 
sequential juxtapositions of amino acid hydrophobicities, using w-point window moving 

20 averages to generate what are commonly known as "hydropathy plots". The largest, 
longest positive variations in these smoothed hydrophobic amplitude graphs across 
sequence-indexed location of membrane proteins are interpreted as the lipophilic, 
hydrophobic transmembrane segments of the membrane protein. The best-studied 
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example of this approach is the finding of seven sequential hydrophobic maxima of 
approximately 25 residues each in the hydropathy plots of bacteriorhodopsin, assumed to 
be the evolutionary prototype of the G-protein gene superfamily of transmembrane 
receptors. This common transmembrane receptor protein motif comprises copolymers of 
5 seven transmembrane domains that snake back and forth across the lipid bilayers of 

membranes, anchored by lipophilic transmembrane ("TM") segments. In this motif, three 
separate extracellular loops ("ELs") are defined by the TMs: the first extracellular loop, 
EL-I, between TM 2 and TM 3 ; the second extracellular loop, EL-II, between TM4 and 
TM 5 ; and the third extracellular loop, EL-III, between TM 6 and TM 7 . 

10 Secondary structures with matching wavenumbers, such as the p-strands of 

interleukin-lp, have been shown to bind together and initiate protein folding in a process 
called the "hydrophobic zipper". We define "wavenumbers" as the inverse spatial 
variational frequencies of a physicochemically transformed series. They are reported 
here in sequential distance units of amino acids. Two long, helical secondary structures 

1 5 with congruent hydrophobic wavenumbers bind to create the central "hydrophobic knot" 
that stabilizes the structure of phospholipase A 2 . Recent studies of the binding of 
extracellular domains of growth hormone receptor by polyclonal antibodies to ovine 
growth hormone have shown that functional binding occurs between the epitope 
sequences and the extracellular segments of the growth hormone transmembrane 

20 receptor. This binding, analogous to that between peptide ligands and their receptors, is 
more related to common helical, loop and/or disordered secondary structures than to 
specific amino acid sequences or their local three-dimensional geometry. 
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Estimates of the relative contributions by the AGn p of each of the twenty amino 
acids to these weak bond-mediated reactions can be approximated as the free energy of 
transfer from aqueous to organic phases of each of the amino acids in a binary solution. 
Values for the free energy of transfer are measured as the relative equilibrium partitions 

5 Ke, = e RT , expressed in kcal/mol, in these aqueous-organic binary solvents. The 

transformation of individual amino acids into their AGh P values enables the conversion of 
polypeptide and protein sequences into real number series available for analyses with 
respect to matches in sequential patterns. These have been predictive of differentially 
selective hydrophobic attraction and aggregation between peptide ligands and relevant 

10 extracellular receptor loops following their search via "snake upon snake" sliding 
diffusion, or "reptation". 

A topologically one-dimensional polypeptide sequence manifests secondary 
structures, which are organized into supersecondary structures and further into tertiary 
structures. For example, spiral rotations of » 3.6 amino acids are the elementary 

1 5 component of a helical barrel comprised of 1 2- 1 6 amino acids. These helical barrels may 
be joined by short loops into four-barrel bundles comprised of 60-70 amino acids, which 
may in turn be part of a protein domain containing several hundred amino acids and 
forming sequentially segregated or alternating barrels, bundles, (5-sheets and coils and 
loops of varying lengths. Therefore, hydrophobic sequences of a range of lengths may 

20 underlie the conformational components of different sizes and complexity that comprise 
the compact intermediate states of proteins. 
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Transformations of polypeptide sequences into AGh P values have been found 
useful in predicting polypeptide chain turns composing secondary structures, such as a- 
helices and p-strands. These predictions have been confirmed by x-ray crystallographic 
studies. Generic ot-helices are « 5.4 angstroms long with 3.6 amino acids per rotation 
5 resulting in * 1 .5 angstrom linear distance per residue. Generic p-strands have 2.1 amino 
acids per turn with * 3.3 angstroms linear distance per residue. 

Sliding window AGh P averages were shown to be able to locate the lipophilic, 
hydrophobic transmembrane segments of membrane proteins, and these results were 
3 confirmed using low- and high-resolution crystallographic studies of bacteriorhodopsin 

10 as a model seven-transmembrane receptor protein. It is generally accepted that 
± representation of polypeptide sequences as a series of amino acid aqueous volumes, 

y partial specific volumes or AGh P , followed by n-block averaging, statistical predilection, 

u hydrophobic moments, Fourier transformation, helical wheel plots or wavelet 

h j transformations can predict the size and locations of secondary and transmembrane 

M- 1 5 structures in soluble and membrane proteins 60-80% of the time. These approaches have 
also been found useful in predicting supersecondary structures, such as the four-helix 
barrels and the supercoiling of a-helical structures about each other in fibrous proteins, 
such as the keratins and myosin tails. However, one drawback of these methods is that 
coexisting sequential variations in hydrophobic free energy wavelengths (mode or 
20 modes) other than that of transmembrane segments are lost in the generation of 

hydropathy plots by smoothing. Moreover, conventional Fourier transformation of the 
protein's hydrophobicities results in poor mode definition, because of end effects and 
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intrinsic multimodality. In addition, these conventional techniques have thus far 
provided no solution of what is called the "inverse problem" - that is, even if the 
conventional methods were able to define one or more given signatory and relevant 
modes, how does one construct a de novo peptide using these modes? The present 
5 invention overcomes the deficiencies of the prior art, and describes successful solutions 
to the inverse problem. 

When the amino acid sequences of neuropeptides and peptide hormones were 
transformed into their individual AGh P values, functionally related peptides demonstrated 
similarities in hydrophobic free energy power spectral mode or modes. Functionally 

10 related peptide family members share the same statistically significant dominant power 
spectral wavelengths (wavenumbers expressed as inverse spatial frequencies), though 
differing in their ordered amino acid content by as much as 60%. The power spectral 
wavelengths are expressed in units of amino acid residues as h{d). For example, 
glucagon, vasoactive intestinal peptide, secretin, oxytomodulin, helodermin and growth 

15 hormone releasing factor, which share several (but not all) physiological actions and 
which have differing relative potencies, share a A(co) = 4.0. The range of peptide 
hydrophobic modes found by the power spectral transformation of amino acid sequences 
as hydrophobic free energies includes the well known h(to) = 3.6 and h(<6) = 2.0 of the a- 
helix and the p-strand, respectively, but many others as well, ranging from the A(oo) = 

20 13.10 amino acid residue of acid fibroblast growth factor to the h(<a) = 2.18 which 

dominates the hydrophobic free energy power spectrum of corticotropin releasing factor. 
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The HIV coat protein manifests a waxing and waning of h(o>) = 7 to 9 (observed 
by sliding a 50-residue windowed Fourier transform along its sequence), which appears 
to be conserved across many of its mutations. Fibroblast growth factor ("FGF") was 
predicted and confirmed to have a regulatory influence on the enzyme ribonuclease A, 
5 with which it was found to share dominant hydrophobic mode. This mode match led to 
experiments that demonstrated an increased half-life of messenger RNA in the presence 
of FGF in a neuroendocrine cell line. 

The specific amino acid sequences of the calcitonins, the peptide hormone family 
that regulates the rate of enzymatic bone catabolism, vary by approximately 60% across 

10 species, but all are dominated by an A(co) = 3.6. The most potent calcitonin (from 

salmon) expresses this mode with a significantly lower hydrophobicity per residue (due 
the presence of a higher number of charged groups) than those of nine other species 
examined. The same h(ti) can be expressed across differing average hydrophobicities of 
the amino acid sequences of peptides and receptors. 

1 5 Using a variety of techniques involving linear decomposition and transformation 

of the AGhp sequences, we have obtained diagnostic graphical patterns of known and 
novel proteins with weak or unknown homology, polyproteins which have multiple 
functional segments following post-translational processing, and discriminable subtypes 
in membrane pore, channel and transporter proteins. These methods, which decompose 

20 AGhp series into their hierarchical levels of organization to yield secondary and 

supersecondary patterns at multiple wavelengths and/or length scales, include a variety of 
wavelet transformations, eigenvalue decomposition of autocovariance matrices and all 
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poles, maximum entropy power spectra. Using AGh P sequences as input, these methods 
elucidated primary and secondary wavenumbers and the sequential order of these 
multiple hydrophobic modes which, when taken together, can contribute to the 
preliminary classification of unknown proteins into families or provide clues to their 
5 function. 

Using these techniques, we have located peptide-receptor mode matches in the 
ELs of seven-transmembrane proteins, in the vicinity of neurotransmitter and 
pharmacological binding domains suggested by studies of point mutations and chimeric 
exchanges. The ligands designed for mode-matched hydrophobic aggregation at these 

10 sites are postulated to have modulatory (e.g. allosteric and/or direct) influences on the 
physiological activities induced by the corresponding membrane protein's native ligands. 
In addition, mode matches were found between the a-estrogen receptor and a known 
peptide antagonist; between a nuclear membrane docking site on a nuclear factor of 
activated T-cells and the known ligand calcineurin; and between the protein chaperonin 

1 5 GroEL and p-lactamase, which is known to be bound by GroEL. 

Eigenfunctions of autocovariance matrices of lagged AGh P sequence data matrices, 
maximum entropy power spectra and wavelet transformations were used as linear 
decompositions to remove the longer AGh P sequence wavelengths of various receptor 
TMs, leaving the shorter wavelength hydrophobic modes for analyses. Matches as 

20 statistical patterns in AGhp modes were found between peptide ligands and their 
membrane receptors, including kappa, mu, delta and orphan opiate receptors, 
corticotropin releasing factor receptor, cholecystokinin receptor, neuropeptide Y receptor, 
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somatostatin receptor, bombesin receptor, and neurotensin receptor. Functionally 
significant mode matches also occur between peptides and non-peptide receptors and 
other proteins. For example, AGh P mode matches, such as those found between the 
dopamine co-localized neuropeptide neurotensin and the D2 dopamine membrane 
5 receptor, D 2 DA, and those found between the gastrointestinal and brain peptide 

cholecystokinin and the dopamine membrane transporter, DAT, predicted the differential 
binding of the pharmacologically active ligands to their respective responsive dopamine 
membrane receptors and, correspondingly, their lack of binding to the opposing, 
pharmacologically unresponsive dopamine membrane receptors. 

10 We have proposed that functional interactions of peptides and biogenic amines 

may occur via selective hydrophobic aggregation of these peptides with mode-matched 
ELs on a target membrane protein. These interactions may result in heterosteric 
modification of the global kinetic conformations of the target membrane protein, and 
thereby produce responses to native or pharmacological ligands, distant from 

1 5 intramembranous ion- or charge-mediated active sites. We have modeled the joint 
actions on a single membrane protein as the shifting of the critical hydrophilic- 
hydrophobic partition between extra- and intramembranous portions of the TMs of 
receptors by peptide-receptor loop hydrophobic weak bond binding. This would facilitate 
(or retard) the first-order phase transition of native ligand induced-receptor membrane 

20 internalization, where low dielectric constant, unscreened ionic and/or charge-mediated 
tight binding most likely occurs. This theory contrasts with another suggesting that 
receptor-mediated interactions between co-localized biogenic amines and neuropeptides, 
such as dopamine and cholecystokinin, result from convergent intramembranous 
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signaling through two receptors, one for each ligand, via the cooperative interactions 
between their membrane receptor proteins which result in G-protein mediated second 
messenger cascades. 

Peptides are known to mediate a variety of physiological responses in many 
organisms, including man. Among these bioactive peptides are the peptide hormones, 
such as glucagon and insulin, which regulate glucose levels in the blood; gastrin and 
secretin, which control digestive processes; and follicle-stimulating hormone (FSH) and 
leuteinizing hormone, which regulate reproductive processes. Other bioactive peptides 
act as growth factors, including somatotropin (growth hormone), erythropoietin, and 
NGF (nerve growth factor). 

Because of the powerful and specific effects of these peptides, they have long 
held great interest as drug candidates. For example, insulin is widely used to combat 
diabetes, and erythropoietin stimulates red blood cell formation. However, peptides have 
numerous drawbacks as potential therapeutics. Peptides are very unstable and sensitive 
to changes in their environments, which can create alterations in their structures and 
reduce or eliminate their physiological effects. Furthermore, peptides are susceptible to 
proteolysis, which complicates the problem of delivery to the desired site in the body and 
limits the available routes of administration. The available routes of administration are 
further limited by the relatively large sizes of many peptides, which make transdermal or 
inhalation administration methods impractical. Because peptides typically interact with 
other peptides or proteins to produce their biological effects, and the in vivo interactions 
between even a simple peptide and another protein are extraordinarily difficult to 
understand, enormous effort is required to determine the interactions between such 
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molecules, or even to predict if such interactions will occur. Finally, relatively few 
bioactive peptides are known, in comparison to the number of potential polypeptide 
targets that mediate biological effects. As a result, there is great interest in finding 
methods to predict sequences of peptides that will interact with a polypeptide/protein 
5 target, and produce a desired physiological response. The present inventors have made 
the revolutionary discovery that peptides, in interaction with solvent-accessible proteins, 
also influence the behavior of proteins (as above) that are not specific peptide receptors. 

The difficulties associated with predicting the structure of peptides that would 
produce a given effect in the body have led to the adoption of various combinatorial 

1 0 approaches. These methods produce large numbers of peptides having randomly 
generated sequences. The peptides are then subjected to various high-throughput 
screening methods to detect those peptides that may warrant further study. However, 
without prior knowledge of a relevant sequence pattern, often called a peptide 
pharmacophore, and without proven methods of pattern-conserving design, finding 

1 5 physiologically active lead compounds in applications involving peptide-protein 

interactions using purely random combinatorial searches is generally a low probability 
event. Depending on the candidate peptide length, the statistical expectations with 
respect to hits in at least micromolar concentrations using high throughput screening of > 
300,000-400,000 component peptide libraries generated by parallel synthesis and 

20 combinatorial strategies, can be less than 2-4 per 100,000 peptides. Detection of these 
candidate peptides requires costly and time-consuming high-throughput methods for both 
peptide synthesis and for screening of the peptides. As a result, there is a great need for a 
method that can produce peptides or peptide-like drugs having a high probability of 
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binding, modulating the activity of, activating or inhibiting a target polypeptide and/or 
protein. 

SUMMARY OF THE INVENTION 
5 The present invention relates to entirely new methods of designing peptides or 

peptide analogue molecules capable of binding to and/or otherwise modulating the 
function of protein targets having known amino acid sequences. The methods employ 
three kinds of templates, derived from analyses of the target protein sequences, in 
addition to relevant distributions of amino acids, for weighted and constrained random 

10 assignments to the templates to produce the peptides. Protein targets suitable for use in 
the present invention include cell membrane receptors, nuclear membrane receptors, 
circulating peptide and non-peptide receptors, membrane and circulating transporters, 
enzymes, chaperonins and chaperonin-like proteins; antibodies, surface proteins of 
infectious agents, and more generally, any protein involved in peptide-protein and/or 

1 5 protein-protein interactions. The peptides are designed to bind to and/or otherwise 

modulate, activate and/or inhibit the function of the target protein. The kinetic influence 
of the algorithmically-designed peptides on target protein function maybe direct, 
competitive, uncompetitive, noncompetitive and/or allosteric in character. The templates 
are derived from at least one of the following: 1) eigenvectors of the autocovariance 

20 matrices of the physicochemically transformed amino acid sequence of the target protein; 
2) wavelet subsequence templates derived from a variety of wavelet transformations of 
the physicochemically transformed amino acid sequence of the target protein; and 3) 
redundant subsequence templates computed from the physicochemically transformed 
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amino acid sequence of the target protein. In the methods of the present invention, the 
constituent amino acids employed in synthesis of the peptide are partitioned into a finite 
number of groups, based on similarities in values of a physicochemical property. 
Thereafter, the amino acids are randomly assigned to the peptide, based on matching the 
physicochemical mode of the template derived from the target protein amino acid 
sequence. Partitioned amino acid distributions for random assignments to the similarly 
partitioned templates may be weighted by, for example, consideration of amino acid 
distribution in a variety of extra- and/or intracellular physiologically relevant pools or 
alternatively, such distributions in regions in the target protein sequence relevant to the 
construction of the templates. The physicochemical transformations of each of the amino 
acids in the target protein sequence may be based on, for example, hydrophobic free 
energy, relative vapor pressure, relative free energy of amino acid transfer into bulk 
phases, aqueous molar volume, aqueous surface area, aqueous cavity surface area, partial 
specific volume, relative charge, relative mass (in daltons), volume, pK a , relative 
diffusivity, relative frictional coefficient, relative chromatographic mobility, relative 
electrophoretic mobility, and/or memberships in categorical amino acid families such as 
polar, uncharged, polar charged, basic-positively charged, acidic-negatively charged and 
sulfur containing. Sequential pattern ("mode") matches between candidate algorithmic 
peptides and their target proteins are designed such that when examined by maximum 
entropy, all poles, power spectral transformations and/or wavelet transformations, they 
yield peaks with wavenumbers that differ by 10% or less of the larger wavenumber value. 
As noted above, wavenumbers are the inverse spatial variational frequencies of a 
physicochemical transformed data series, expressed in sequential distance units of amino 
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acids. These peptides are then selected for physiological testing on the target protein 
system. The peptide design methods and an associated mechanistic rationale are 
illustrated for the methods of the present invention, using an eigenvector template derived 
from the hydrophobic free energy-transformed sequence of several different receptors 

5 and random assignment of amino acids to the eigenvector templates based on probability- 
weighted amino acid pool distributions. The peptides generated in this manner 
demonstrate physiological activity in receptor-transfected cell systems, as shown by 
direct action and/or pretreatment potentiation or inhibition of extracellular acidification 
rates. In addition, peptides generated by the methods of the present invention also bound 

10 to and otherwise interact with and alter the activities of the seven-transmembrane 

cholinergic Ml receptor ("muscarinic Ml receptor") and the nerve growth factor (NGF) 
receptor, which has one transmembrane segment. As another example of the range of 
applicability of these methods, hydrophobic free energy mode matches between the 
peptide fibroblastic growth factor and ribonuclease successfully predicted their functional 

1 5 interaction in neuroendocrine cell culture. These results illustrate the broad applicability 
of the methods of the present invention to the design of peptides for binding to or 
otherwise modulating a wide variety of different kinds of target polypeptides and 
proteins. 

One of the three mode-matched peptide design methods of the invention involves 
20 the construction of such peptides using random assignment of peptide constituents, such 
as amino acids, as dictated by an eigenvector template containing polypeptide-matching 
physicochemical property binding/modulating modes. This method is herein exemplified 
by one of many possible physicochemical properties usable in the method, namely, 
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hydrophobic free energy. The template eigenvector is obtained by linear decomposition 
of an autocovariance matrix formed by transformation of the polypeptide's amino acid 
sequence into a physicochemical sequence, in this case a hydrophobic free energy data 
series. The leading eigenvalue-associated eigenvectors are convolved with the original 

5 hydrophobic free energy data series to construct eigenfunctions. These eigenfunctions 
may then be further analyzed using wavelet transformations and all poles, maximum 
entropy power spectral transformations. The wavelet transformations may be discrete or 
continuous, and further may be one-dimensional wavelet packets or multiple convolved 
wavelet transformations. This approach yields clean representations of the polypeptide 

10 hydrophobic free energy modes as leading and secondary eigenfunctions. Most of the 
information found in the secondary eigenfunctions would be lost in the conventional 
smoothing of hydropathy plots, or contaminated by end effects and multimodality in 
conventional Fourier transformations. The eigenvectors associated with these 
eigenfunctions are used as templates for the formation of mode-matched peptides that can 

1 5 be tested for their ability to bind to or otherwise modulate the receptor. A mode match is 
attained when the maximum entropy power spectral or wavelet transformations of the 
polypeptide and the peptide or peptide-like molecule yield wavenumbers that differ by 
10% or less of the larger wavenumber value. The amino acids intended for use in 
producing the candidate peptide are grouped into a number of groups, based on their 

20 assigned values of a physicochemical property (e.g. hydrophobic free energy). The 
eigenvector associated with the eigenfunction (or, alternately, the eigenvectors-based 
vector) is graphed, where the x-axis shows ordered position of the eigenvector and the Y- 
axis shows the numerical values of the physicochemical property. The y-axis is 
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partitioned into an equal number of groups as intervals of the y-axis (e.g., four equal 
intervals), converting the eigenvector (or eigenvectors-based vector) into an eigenvector 
template. Amino acids corresponding to the value of the physicochemical property on 
the y-axis of the eigenvector template are randomly assigned to positions in the template, 
5 forming peptides or peptide-like molecules. The amino acid assignments may also be 
weighted or otherwise altered in accordance with a specific amino acid pool distribution 
or in accordance with known effects of substitutions of individual amino acids or amino 
acid segments, if desired. 

The second method involves the construction of mode-matched peptides through 

10 the generation of wavelet subsequence templates derived from a variety of wavelet 

transformations of the physicochemically-transformed amino acid sequence of the target 
protein. The wavelet transformation method is particularly well suited for the study of 
localized coherent structures that appear across a target protein sequence, such as the 
patterns of alternating helices, loops and strands that make up larger supersecondary 

1 5 structures, such as helical barrels and sheets. A number of mother wavelet families are 
available for use in wavelet transformations. 

The third method produces redundant target polypeptide or protein subsequence 
templates from the physicochemically-transformed amino acid sequence of the target 
polypeptide or protein. Redundant subsequence templates are prepared by converting the 

20 amino acid sequence of the target polypeptide or protein into a template through 

symbolic representations of each amino acid, e.g., one-letter amino acid codes or, more 
preferably, values representing each amino acid's membership in a particular 
physicochemical property grouping. The transformed target polypeptide or protein 



20 



PATENT 

EXPRESS MAIL NO. EL557885546US 
sequence is then scanned to find all possible redundant nonoverlapping subsequences. 
The redundant subsequences detected are used as templates to create mode-matching 
peptides. 

It is therefore an object of the present invention to provide a method for 
5 synthesizing a peptide or a peptide-like molecule based on matching a physicochemical 
mode of a target polypeptide or protein to the same physicochemical mode of the peptide 
or peptide-like molecule, comprising the steps of assigning a numerical value of an 
orderable physicochemical property to each member of a set of peptide constituents 
which includes all the members of the set of naturally-occurring amino acids, arranging 

1 0 the peptide constituents in order of the numerical values of an orderable physicochemical 
property, partitioning the set of peptide constituents into a plurality of peptide constituent 
groups, whereby each of the peptide constituent groups contains at least one member of 
the set of peptide constituents, each peptide constituent group encompasses a range of the 
numerical values, each member of the set of peptide constituent belongs to only one 

1 5 peptide constituent group, creating a polypeptide physicochemical data series by 

replacing each amino acid in an amino acid sequence of the target polypeptide or protein 
with the numerical value of the orderable physicochemical property corresponding to 
each amino acid in the amino acid sequence, calculating one or more polypeptide 
eigenvalues and a corresponding polypeptide eigenvector associated with each of the 

20 polypeptide eigenvalues by linear decomposition of an autocovariance matrix formed 
from a sequentially lagged data matrix of the polypeptide physicochemical data series, 
ordering the polypeptide eigenvalues and the corresponding polypeptide eigenvectors 
from largest to smallest, selecting one or more of the polypeptide eigenvectors, 
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transforming the selected polypeptide eigenvectors into an eigenvector template, forming 
a graph of the eigenvector template, wherein the numerical values of the physicochemical 
property are graphed along the y-axis of the graph and ordered position in the eigenvector 
template is graphed along the x-axis of the graph, partitioning the graph along the y-axis 
5 according to the ranges of the numerical values of the physicochemical property defining 
the peptide constituent groups to form a plurality of y-axis ranges, assigning a member of 
the peptide constituent group to each position in the peptide or peptide-like molecule by 
using the graph as a template, wherein at each ordered position in the eigenvector 
template along the x-axis of the graph, the member of the peptide constituent group 

1 0 assigned to the ordered position has a value of the orderable physicochemical property 
that is within the y-axis range of the ordered point, and synthesizing the peptide or 
peptide-like molecule. 

It is another object of the present invention to provide a method for matching a 
physicochemical mode of a peptide or a peptide-like molecule to the same 

1 5 physicochemical mode of a target polypeptide or protein to determine if the peptide will 
bind to and/or otherwise modulate the target polypeptide or protein, comprising the steps 
of assigning a numerical value of an orderable physicochemical property to each member 
of a set of peptide constituents which includes all the members of the set of naturally- 
occurring amino acids, arranging the peptide constituents in order of the numerical values 

20 of the orderable physicochemical property, partitioning the set of peptide constituents 
into a plurality of peptide constituent groups, whereby each of the peptide constituent 
groups contains at least one member of the set of peptide constituents, each peptide 
constituent group encompasses a range of the numerical values, each member of the set 
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of peptide constituents belongs to only one peptide constituent group, creating a 
polypeptide physicochemical data series by replacing each amino acid in an amino acid 
sequence of the target polypeptide or protein with the numerical value of the orderable 
physicochemical property corresponding to each amino acid in the amino acid sequence, 
5 calculating one or more polypeptide eigenvalues and a corresponding polypeptide 

eigenvector associated with each of the polypeptide eigenvalues by linear decomposition 
of an autocovariance matrix formed from a sequentially lagged data matrix of the 
polypeptide physicochemical data series, ordering the polypeptide eigenvalues and the 
corresponding polypeptide eigenvectors from largest to smallest, transforming the 

1 0 polypeptide physicochemical data series into one or more polypeptide eigenfunctions, 
using the ordered polypeptide eigenvectors as multiplicative weights, transforming the 
polypeptide eigenfunctions into dominant wavenumbers, using all poles maximum 
entropy power spectra, to produce polypeptide spectral power ,peaks, identifying the 
polypeptide power spectral peaks, creating a peptide physicochemical data series by 

1 5 replacing each peptide constituent in a peptide sequence of the peptide or a peptide-like 
molecule with the numerical value of the orderable physicochemical property 
corresponding to the peptide constituent in the peptide sequence, calculating one or more 
peptide eigenvalues and a corresponding peptide eigenvector associated with each of the 
peptide eigenvalues by linear decomposition of an autocovariance matrix formed from 

20 the peptide physicochemical data series, ordering the peptide eigenvalues and the 
corresponding eigenvectors from largest to smallest, transforming the peptide 
physicochemical data series into one or more peptide eigenfunctions, using the ordered 
peptide eigenvectors as multiplicative weights, transforming the peptide eigenfunctions 
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into dominant wavenumbers, using all poles maximum entropy power spectra, to produce 
peptide spectral power peaks, identifying the peptide power spectral peaks, and 
comparing the polypeptide spectral power peaks to the peptide spectral power peaks to 
determine if the polypeptide spectral power peaks match the peptide spectral power 
5 peaks, wherein a match between the polypeptide spectral power peaks and the peptide 
spectral power peaks indicates the peptide or peptide-like molecule may bind to and/or 
otherwise modulate the target polypeptide or protein. 

It is another object of the present invention to provide a method for matching a 
p peptide or a peptide-like molecule to a target polypeptide or protein to determine if the 

% i 1 0 peptide will bind to and/or otherwise modulate the target polypeptide or protein, 
% i comprising the steps of assigning a numerical value of an orderable physicochemical 

if; property to each member of a set of peptide constituents, the set of peptide constituents 

JU including all the members of the set of naturally-occurring amino acids, arranging the 

nj peptide constituents in order of the numerical values of the orderable physicochemical 

O 1 5 property, partitioning the set of peptide constituents into a plurality of peptide constituent 
groups, whereby each of the peptide constituent groups contains at least one member of 
the set of peptide constituents, each peptide constituent group encompasses a range of the 
numerical values, each member of the set of peptide constituents belongs to only one 
peptide constituent group, creating a polypeptide physicochemical data series by 
20 replacing each amino acid in an amino acid sequence of the target polypeptide or protein 
with the numerical value corresponding to the amino acid in the amino acid sequence, 
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decomposing the polypeptide physicochemical data series into translated and scaled 

' i — b 

version of a mother wavelet, w, as V^fab) ={\l4a) J#0'M )di 

o a 

wherein w denotes the chosen mother wavelet function, separating Vf^{a,b) into 
polypeptide modulus and polypeptide phase parts, graphing the polypeptide phase parts 

5 on a polypeptide phase graph, wherein the x-axis of the polypeptide phase graph indexes 
sequence position and the y-axis of the polypeptide phase graph is numbered in units of 
one of dilate divisions (dd) and wavelet wavelengths (gj), graphing the polypeptide 
modulus parts on a polypeptide modulus graph, wherein the x-axis of the polypeptide 
modulus graph indexes sequence position and the y-axis of the polypeptide modulus 

1 0 graph is numbered in units of one of dilate divisions (dd) and wavelet wavelengths (to), 
identifying a plurality of polypeptide maximal phase amplitudes and a plurality of 
polypeptide moduli in the polypeptide phase graph and the polypeptide modulus graph, 
respectively, creating a peptide physicochemical data series by replacing each peptide 
constituent in a peptide sequence of the peptide or a peptide-like molecule with the 

1 5 numerical value of the orderable physicochemical property corresponding to each the 
peptide constituent in the peptide sequence, decomposing the peptide physicochemical 
data series into translated and scaled version of a mother wavelet, w, as V^(a,b) = 

(l/Ja'))H(i)w(—)di 
o a 

wherein w denotes the chosen mother wavelet function, separating J^(a,b) into peptide 
20 modulus and peptide phase parts, graphing the peptide phase parts on a peptide phase 

graph, wherein the x-axis of the peptide phase graph indexes sequence position and the y- 
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axis of the peptide phase graph is numbered in units of one of relative dilation (dd) and 
wavelet wavelengths (to), graphing the peptide modulus parts on a peptide modulus 
graph, wherein the x-axis of the peptide modulus graph indexes sequence position and the 
y-axis of the peptide modulus graph is numbered in units of one of dilate divisions (dd) 

5 and wavelet wavelengths (to), identifying a plurality of peptide maximal phase 

amplitudes and a plurality of peptide moduli in the peptide phase graph and the peptide 
modulus graph, respectively, comparing the plurality of polypeptide maximal phase 
amplitudes in the polypeptide phase graph to the plurality of peptide maximal phase 
amplitudes in the peptide phase graph to determine if the plurality of polypeptide 

10 maximal phase amplitudes match the plurality of peptide maximal phase amplitudes, 
comparing the plurality of polypeptide moduli in the polypeptide modulus graph to the 
plurality of peptide moduli in the peptide modulus graph to determine if the plurality of 
polypeptide moduli match the plurality of peptide moduli, wherein a match between the 
plurality of polypeptide maximal phase amplitudes and the plurality of peptide maximal 

15 phase amplitudes, and a match between the plurality of polypeptide moduli and the 
plurality of peptide moduli, indicates the peptide or peptide-like molecule may bind to 
and/or otherwise modulate the polypeptide. 

It is another object of the present invention to provide a method for matching a 
peptide or a peptide-like molecule to a target polypeptide or protein to determine if the 

20 peptide will bind to and/or otherwise modulate the target polypeptide or protein, 

comprising the steps of assigning a numerical value of an orderable physicochemical 
property to each member of a set of peptide constituents, the set of peptide constituents 
including all the members of the set of naturally-occurring amino acids, arranging the 
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peptide constituents in order of the numerical values of the orderable physicochemical 
property, partitioning the set of peptide constituents into a plurality of peptide constituent 
groups, whereby each of the peptide constituent groups contains at least one member of 
the set of peptide constituents, each group encompasses a range of the numerical values, 
each member of the set of peptide constituents belongs to only one peptide constituent 
group, creating a polypeptide physicochemical data series by replacing each amino acid 
in an amino acid sequence of the target polypeptide or protein with the numerical value 
corresponding to the amino acid in the amino acid sequence, decomposing the 
polypeptide physicochemical data series with a family of functions 
W j „i ! (x)=2~ jn W tt (2~ J x-k) , which when j,n are positive integers and k has an integer 
value, are organized in one or more tree structures, each of the tree structures being 
composed of a plurality of nodes, each of the nodes being in the form of: 

A 

wherein Wj, „,*(x) is computed for a mother wavelet function, computing and frequency 
ordering best level and best tree representations of a physicochemical polypeptide series 
based on Stein's Unbiased Risk Estimate (SURE) and Shannon entropy criteria, graphing 
the best level representation on a polypeptide best level graph, wherein the x-axis of the 
polypeptide best level graph indexes sequence position and the y-axis of the polypeptide 
best level graph is numbered in units of wavelet wavelengths, tn , graphing the best tree 
representation on a polypeptide best tree graph, wherein the x-axis of the polypeptide best 
tree graph indexes sequence position and the y-axis of the polypeptide best tree graph is 
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numbered in units of one of relative dilation (dd) and wavelet wavelengths, xn , 
identifying a plurality of polypeptide maximal coefficient amplitudes, each of the 
plurality of polypeptide maximal coefficient amplitudes being derived from the 
polypeptide best level graph and the polypeptide best tree graph, creating a peptide 
physicochemical data series by replacing each peptide constituent in a peptide sequence 
of the peptide or a peptide-like molecule with the numerical value of the orderable 
physicochemical property corresponding to the peptide constituent in the peptide 
sequence, decomposing the peptide physicochemical data series with the family of 
functions W J ,„, k ix)=2' jn W n (2' J x-k), which when/,H are positive integers and & has an 
integer value, are organized in one or more tree structures, each of the tree structures 
being comprised of a plurality of nodes, each of the nodes being in the form of 

A 

Wj+l,2n Wj+j.2n+l 

wherein ^,„.*( x ) is computed for a mother wavelet function, computing and frequency 
ordering best level and best tree representations of a physicochemical peptide series 
based on SURE and Shannon entropy criteria, graphing the best level representation on a 
peptide best level graph, wherein the x-axis of the peptide best level graph indexes 
sequence position and the y-axis of the peptide best level graph is numbered in units of 
wavelet wavelengths, m , graphing the best tree representation on a peptide best tree 
graph, wherein the x-axis of the peptide best tree graph indexes sequence position and the 
y-axis of the peptide best tree graph is numbered in units of one of relative dilation (dd) 
and wavelet wavelengths, m , identifying a plurality of peptide maximal coefficient 
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amplitudes, each of the plurality of peptide maximal coefficient amplitudes being derived 
from the peptide best level graph and the peptide best tree graph, comparing the plurality 
of polypeptide maximal coefficient amplitudes to the plurality of peptide maximal 
coefficient amplitudes to determine if the plurality of polypeptide maximal coefficient 
5 amplitudes match the plurality of peptide maximal coefficient amplitudes, wherein a 
match between the plurality of polypeptide maximal coefficient amplitudes and the 
plurality of peptide maximal coefficient amplitudes indicates the peptide or peptide-like 
molecule may bind to and/or otherwise modulate the target polypeptide or protein. 

It is another object to provide a method for modifying a non-peptide-responsive 

1 0 target polypeptide or protein to bind to and/or otherwise modulate a peptide or peptide- 
like molecule by modifying the sequence of the non-peptide-responsive target 
polypeptide or protein to match a physicochemical mode of the peptide or peptide-like 
molecule, comprising the steps of assigning a numerical value of an orderable 
physicochemical property to each member of a set of polypeptide constituents, the set of 

1 5 peptide constituents including all the members of the set of naturally-occurring amino 

acids, arranging the peptide constituents in order of the numerical values of the orderable 
physicochemical property, partitioning the set of peptide constituents into a plurality of 
peptide constituent groups, whereby each of the peptide constituent groups contains at 
least one member of the set of peptide constituents, each group encompasses a range of 

20 the numerical values, each member of the set of peptide constituents belongs to only one 
peptide constituent group, creating a polypeptide physicochemical data series by 
replacing each amino acid in an amino acid sequence of the non-peptide-responsive target 
polypeptide or protein with the numerical value corresponding to the amino acid in the 
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amino acid sequence, calculating one or more polypeptide eigenvalues and a 
corresponding polypeptide eigenvector associated with each of the polypeptide 
eigenvalues by linear decomposition of an autocovariance matrix formed from the 
polypeptide physicochemical data series, ordering the polypeptide eigenvalues and the 
5 corresponding polypeptide eigenvectors from largest to smallest, transforming the 
polypeptide physicochemical data series into polypeptide eigenfunctions, using the 
ordered polypeptide eigenvectors as multiplicative weights, transforming the polypeptide 
eigenfunctions into dominant wavenumbers, using all poles maximum entropy power 
spectra to produce polypeptide spectral power peaks, identifying the polypeptide power 

10 spectral peaks, creating a peptide physicochemical data series by replacing each peptide 
constituent in a peptide sequence of the peptide or peptide-like molecule with a numerical 
value of the orderable physicochemical property corresponding to the peptide or peptide- 
like molecule constituent in the peptide sequence, calculating one or more peptide 
eigenvalues and a corresponding peptide eigenvector associated with each of the peptide 

1 5 eigenvalues by linear decomposition of an autocovariance matrix formed from the 
peptide physicochemical data series, ordering the peptide eigenvalues and the 
corresponding peptide eigenvectors from largest to smallest, transforming the peptide 
physicochemical data series into peptide eigenfunctions, using the peptide eigenvectors 
as multiplicative weights, transforming the peptide eigenfunctions into dominant 

20 wavenumbers, using all poles maximum entropy power spectra, to produce peptide 
spectral power peaks, identifying the peptide power spectral peaks, comparing the 
polypeptide spectral power peaks to the peptide spectral power peaks to determine if the 
polypeptide spectral power peaks match the peptide spectral power peaks, wherein a 
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match between the polypeptide spectral power peaks and the peptide spectral power 
peaks indicates the peptide or peptide-like molecule may bind to and/or otherwise 
modulate the non-peptide-responsive target polypeptide or protein, and if the polypeptide 
spectral power peaks do not match the peptide spectral power peaks, modifying the 
5 amino acid sequence of the non-peptide-responsive target polypeptide or protein to form 
a match between the polypeptide spectral power peaks and the peptide spectral power 
peaks. 

It is a further object to provide a method for modifying a non-peptide-responsive 
5 target polypeptide or protein to bind to and/or otherwise modulate a peptide or peptide- 

^ 10 like molecule by modifying the sequence of the non-peptide-binding/modulating target 
5 ; polypeptide to match a physicochemical mode of the peptide or peptide-like molecule, 

% comprising the steps of assigning a numerical value of an orderable physicochemical 

O property to each member of a set of peptide constituents, the set of peptide constituents 

I U including all the members of the set of naturally-occurring amino acids, arranging the 

3 1 5 peptide constituents in order of the numerical values of the orderable physicochemical 

property, partitioning the set of peptide constituents into a plurality of peptide constituent 
groups, whereby each of the peptide constituent groups contains one or more members of 
the set of peptide constituents, each group encompasses a range of said numerical values, 
each member of the set of peptide constituents belongs to only one peptide constituent 
20 group, creating a polypeptide physicochemical data series by replacing each amino acid 
in an amino acid sequence of the non-peptide-binding and/or modulating target 
polypeptide or protein with a numerical value corresponding to each the amino acid in the 
amino acid sequence, decomposing the polypeptide physicochemical data series into 
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translated and scaled version of a mother wavelet, w, as 

(l/Va)J#(/M— ¥i 

o a 

wherein w denotes the chosen mother wavelet function, separating ^(a.b) into 
polypeptide modulus and polypeptide phase parts, graphing the polypeptide phase parts 

5 on a polypeptide phase graph, wherein the x-axis of the polypeptide phase graph indexes 
sequence position and the y-axis of the polypeptide phase graph is numbered in units of 
one of relative dilation (dd) and wavelet wavelengths (w), graphing the polypeptide 
modulus parts on a polypeptide modulus graph, wherein the x-axis of the polypeptide 
modulus graph indexes sequence position and the y-axis of the polypeptide modulus 

1 0 graph is numbered in units of one of relative dilation (dd) and wavelet wavelengths (tu), 
identifying a plurality of polypeptide maximal phase amplitudes and a plurality of 
polypeptide moduli in the polypeptide phase graph and the polypeptide modulus graph, 
respectively, creating a peptide physicochemical data series by replacing each peptide 
constituent in a peptide sequence of a peptide or a peptide-like molecule with the 

1 5 numerical value corresponding to each peptide constituent in the peptide sequence, 

decomposing the peptide physicochemical data series into translated and scaled version 

1 i — b 

of a mother wavelet, w, as ^(a.b) = (l/y[a) J#(i'M )di 

o a 

wherein w denotes the chosen mother wavelet function, separating ^(a.b) into peptide 
modulus and peptide phase parts, graphing the peptide phase parts on a peptide phase 
20 graph, wherein the x-axis of the peptide phase graph indexes sequence position and the y- 
axis of the peptide phase graph is numbered in units of one of relative dilation (dd) and 
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wavelet wavelengths (w), graphing the peptide modulus parts on a peptide modulus 
graph, wherein the x-axis of the peptide modulus graph indexes sequence position and the 
y-axis of the peptide modulus graph is numbered in units of one of relative dilation (dd) 
and wavelet wavelengths (tu), identifying a plurality of peptide maximal phase 
5 amplitudes and a plurality of peptide moduli in each of the peptide phase graph and the 
peptide modulus graph, respectively, comparing the plurality of polypeptide maximal 
phase amplitudes in the polypeptide phase graph to the plurality of peptide maximal 
phase amplitudes in the peptide phase graph respectively to determine if the plurality of 
polypeptide maximal phase amplitudes match the plurality of peptide maximal phase 

10 amplitudes, comparing the plurality of polypeptide moduli in the polypeptide modulus 
graph to the plurality of peptide moduli in the peptide modulus graph to determine if the 
plurality of polypeptide moduli match the plurality of peptide moduli, wherein a match 
between the plurality of polypeptide maximal phase amplitudes and the plurality of 
peptide maximal phase amplitudes, and a match between the plurality of polypeptide 

1 5 moduli and the plurality of peptide moduli indicates the peptide or peptide-like molecule 
may bind to and/or otherwise modulate the non-peptide-binding and/or modulating target 
polypeptide or protein, and if the plurality of polypeptide maximal phase amplitudes do 
not match the plurality of peptide maximal phase amplitudes, or if the plurality of 
polypeptide moduli do not match the plurality of peptide moduli, modifying the amino 

20 acid sequence of the non-peptide-binding and/or modulating target polypeptide or protein 
to form a match between the plurality of polypeptide maximal phase amplitudes and the 
plurality of peptide maximal phase amplitudes, and between the polypeptide moduli and 
the peptide moduli. 
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It is a further object to provide a method for modifying a non-peptide-responsive 
target polypeptide or protein to bind to and/or otherwise modulate a peptide or peptide- 
like molecule by modifying the sequence of the non-peptide-responsive target 
polypeptide or protein to match a physicochemical mode of the peptide or peptide-like 
5 molecule, comprising the steps of assigning a numerical value of an orderable 

physicochemical property to each member of a set of peptide constituents, the set of 
peptide constituents including all the members of the set of naturally-occurring amino 
acids, arranging the peptide constituents in order of the numerical values of the orderable 
physicochemical property, partitioning the set of peptide constituents into a plurality of 

10 peptide constituent groups, whereby each of the peptide constituent groups contains one 
or more members of the set of peptide constituents, each group encompassing a range of 
said numerical values, each member of the set of peptide constituents belongs to only one 
peptide constituent group, creating a polypeptide physicochemical data series by 
replacing each amino acid in an amino acid sequence of the non-peptide-binding and/or 

1 5 modulating target polypeptide or protein with the numerical value of the orderable 

physicochemical property corresponding to the amino acid in the amino acid sequence, 
decomposing the polypeptide physicochemical data series with a family of functions 
Wj nJ ix)= 2~ jil W n (2~ J x-k), which when j,n are positive integers and k has an integer 
value, are organized in one or more tree structures, each of the tree structures being 

20 comprised of a plurality of nodes, each of the nodes being in the form of: 

A 

Wj+ijn Wj+1.2n+l 
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wherein the 0y,„.*(x) is computed for a mother wavelet function, computing and 
frequency ordering best level and best tree representations of the physicochemical 
polypeptide series based on SURE and Shannon entropy criteria, graphing the best level 
representation on a polypeptide best level graph, wherein the x-axis of the polypeptide 

5 best level graph indexes sequence position and the y-axis of the polypeptide best level 
graph is numbered in units of wavelet wavelengths, vj , graphing the best tree 
representation on a polypeptide best tree graph, wherein the x-axis of the polypeptide best 
tree graph indexes sequence position and the y-axis of the polypeptide best tree graph is 
numbered in units of one of relative dilation (dd) and wavelet wavelengths, cr , 

10 identifying a plurality of polypeptide maximal coefficient amplitudes, each of the 
plurality of polypeptide maximal coefficient amplitudes being derived from the 
polypeptide best level and best tree graphs, decomposing the peptide physicochemical 
data series with a family of functions Wj,„ ik {x)=2- j,2 W n (2- J x-k) , which wheny,n are 
positive integers and £ has an integer value, are organized in one or more tree structures, 

1 5 each of the tree structures being comprised of a plurality of nodes, each of the nodes 
being in the form of: 

A 

Wj+i,2n Wj+i_2n+l 

wherein the Wj^x) is computed a mother wavelet function, computing and frequency 
20 ordering best level and best tree representations of the physicochemical peptide series 
based on SURE and Shannon entropy criteria, graphing the best level representation on a 
peptide best level graph, wherein the x-axis of the peptide best level graph indexes 
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sequence position and the y-axis of the peptide best level graph is numbered in units of 
wavelet wavelengths, m , graphing the best tree representation on a peptide best tree 
graph, wherein the x-axis of the peptide best tree graph indexes sequence position and the 
y-axis of the best tree graph is numbered in units of one of relative dilation (dd) and 
5 wavelet wavelengths, m , identifying a plurality of peptide maximal coefficient 

amplitudes, each of the plurality of peptide maximal coefficient amplitudes being derived 
from the peptide best level and best tree graphs, comparing the plurality of polypeptide 
moduli in the polypeptide modulus graph to the plurality of peptide moduli in the peptide 
modulus graph to determine if the plurality of polypeptide moduli match the plurality of 

10 peptide moduli, wherein a match between the plurality of polypeptide maximal phase 
amplitudes and the plurality of peptide maximal phase amplitudes, and a match between 
the plurality of polypeptide moduli and the plurality of peptide moduli indicates the 
peptide or peptide-like molecule may bind to and/or otherwise modulate the non-peptide- 
binding and/or modulating target polypeptide or protein, and if the plurality of 

1 5 polypeptide maximal phase amplitudes do not match the plurality of peptide maximal 
phase amplitudes, or if the plurality of polypeptide moduli do not match the plurality of 
peptide moduli, modifying the amino acid sequence of the non-peptide-binding and/or 
modulating target polypeptide or protein to form a match between the plurality of 
polypeptide maximal phase amplitudes and the plurality of peptide maximal phase 

20 amplitudes, and between the polypeptide moduli and the peptide moduli. 

The present invention also provides a method of detecting a cancerous cell or 
tissue, comprising contacting all or a portion of the cancerous cell or tissue with an 
effective amount of a peptide or peptide-like molecule having a physicochemical mode 
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that matches a physicochemical mode of a target polypeptide or protein found on the 
cancerous cell or tissue. 

The present invention also provides a method of detecting a tumor in a patient, 
comprising administering to the patient an effective amount of a peptide or peptide-like 
molecule having a physicochemical mode that matches a physicochemical mode of a 
polypeptide or protein found on the tumor, and detecting binding and/or modulating of 
the peptide or peptide-like molecule to the polypeptide or protein. 

The present invention also provides a pharmaceutical composition for treatment 
of a tumor, comprising a peptide or peptide-like molecule having a physicochemical 
mode that matches a physicochemical mode of a polypeptide or protein found on the 
tumor, and a pharmaceutically acceptable carrier. 

The present invention also provides a diagnostic kit for use in detecting a 
polypeptide or protein, comprising a container having a peptide or peptide-like molecule, 
the peptide or peptide-like molecule having a physicochemical mode that matches a 
physicochemical mode of the polypeptide or protein. 

The present invention also provides a method for screening for a disease 
condition, comprising contacting a sample obtained from a patient with an effective 
amount of a peptide or peptide-like molecule having a physicochemical mode that 
matches a physicochemical mode of a polypeptide or protein found in the sample, 
wherein the presence, absence or abnormality in the polypeptide or protein is diagnostic 
of the presence of the disease condition. 

The present invention also provides a method for screening a member selected 
from the group consisting of water, food, and soil for the presence of a contaminant, 
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comprising contacting the member with a peptide or peptide-like molecule having a 
physicochemical mode that matches a physicochemical mode of a polypeptide or protein 
found in the member, wherein the presence, absence, or abnormality in the polypeptide or 
protein is diagnostic of the presence of the contaminant. 

The present invention also provides a method for treating a disease condition, 
comprising administering to a patient in need of such treatment a peptide or peptide-like 
molecule having a physicochemical mode that matches a physicochemical mode of a 
polypeptide or protein found in the sample, wherein the peptide or peptide-like molecule 
is capable of effecting a direct action and/or modulation of an activity of the polypeptide 
or protein, and the direct action and/or modulation effected by the peptide or peptide-like 
molecule is associated with a change in the disease condition. 

The present invention also provides a method for detecting an interaction between 
a peptide and a target polypeptide or protein, comprising incubating a peptide prepared 
by at least one of the methods of the present invention with the target polypeptide or 
protein under conditions that promote the interaction of the peptide with the target 
polypeptide or protein, and detecting the interaction of the peptide with the target 

polypeptide or protein. 

The present invention also provides a pharmaceutical composition for treatment 
of a disease condition, comprising a peptide or peptide-like molecule having a 
physicochemical mode that matches a physicochemical mode of a polypeptide or protein 
found in the sample, the peptide or peptide-like molecule being capable of effecting a 
direct action and/or modulation of an activity of the polypeptide or protein, and the direcl 
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action and/or modulation effected by the peptide or peptide-like molecule is associated 
with a change in the disease condition, and a pharmaceutically acceptable carrier. 

The above and other objects, features and advantages of the present invention will 
become apparent from the following description read in conjunction with the 
accompanying drawings. 

TTRTEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a flowchart which summarizes the methods of the present invention. 

Figure 2 A (left) is a graph of the hydrophobic free energy series, H t , of the human 
D 2 DA receptor and (right) its broad band, multimodal all poles, maximum entropy power 
spectral transformation iS(co). 

Figure 2B (left) is a graph of the human D 2 DA receptor's dominant eigenfunction, 
Ti, demonstrating the « 7 peaks characteristic of the leading receptor eigenfunction of 
members of the seven-transmembrane receptor superfamily and (right) the associated 
long wavelength peak (> 50 residues) in the S(g>). 

Figure 2C (left) is a graph of the human D 2 DA receptor's secondary 
eigenfunction, ¥2, and (right) its associated peaks in the 5(co) at wavelengths of 8.12 and 
2.61 residues. 

Figure 2D (left) is a graph of the human D 2 DA receptor's secondary eigenvector, 
X 2 , used in the design of new peptides, and (right) its associated peaks in the 5(o) at 
wavelengths of 8.16 and 2.67 residues. 
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Figure 3A is a graph of the wavelet subspace transformation of the H t of the 
D2DA receptor, wherein tn ^Jidd) = 2.3 residues. Sequence position is graphed along 
the x-axis and phase amplitude along the y-axis. 

Figure 3B is a graph of the wavelet subspace transformation of the Hi of the 
5 D2DA receptor, wherein m =ftdd) =8.1 residues. Sequence position is graphed along 
the x-axis and phase amplitude along the y-axis. 

Figure 4A is a graph showing the effects of the SHQR peptide (SEQ ID NO: 1) on 
the EAR responses of the human D2DA-transfected mouse LtK cell system to dopamine 
infusion. DA = control with dopamine alone. 
10 Figure 4B is a graph showing the effects of the THQA (SEQ ID NO:2) peptide on 

the EAR responses of the human D2DA-transfected mouse LtK cell system to dopamine 
infusion. DA = control with dopamine alone. 

Figure 4C is a graph showing the effects of the SHQR (SEQ ID NO:l) peptide on 
the EAR responses of the human D2DA-transfected mouse CHO cell system to dopamine 
15 infusion. DA = control with dopamine alone. 

Figure 4D is a graph showing the effects of the THQA (SEQ ID NO:2) peptide on 
the EAR responses of the human D2DA-transfected mouse CHO cell system to dopamine 
infusion. DA = control with dopamine alone. 

Figure 5 A is a graph showing the effects of the E. . .PL (SEQ ID NO:3) peptide on 
20 the EAR responses of the human D2DA-transfected mouse LtK cell system to dopamine 
infusion. DA = control with dopamine alone. 
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Figure 5B is a graph showing the effects of the E. . .PY (SEQ ID NO:4) peptide on 
the EAR responses of the human D 2 DA-transfected mouse LtK cell system to dopamine 
infusion. DA = control with dopamine alone. 

Figure 5C is a graph showing the effects of the E. . .PL (SEQ ID NO:3) peptide on 
the EAR responses of the human D 2 DA-transfected mouse CHO cell system to dopamine 
infusion. DA = control with dopamine alone. 

Figure 5D is a graph showing the effects of the E. . .PY peptide (SEQ ID NO:4) on 
the EAR responses of the human D 2 DA-transfected mouse CHO cell system to dopamine 
infusion. DA = control with dopamine alone. 

Figure 6A is a graph showing the effects of the Ml receptor-derived peptide ITFT 
(SEQ ID NO:9) on the EAR responses of the human Ml receptor-transfected CHO cell 
system to carbachol infusion, left, control with carbachol alone, right, carbachol plus 
ITFT peptide. 

Figure 6B is a graph showing the effects of the Ml receptor-derived peptide 
FSFQ (SEQ ID NO:7) on the EAR responses of the human Ml receptor-transfected CHO 
cell system to carbachol infusion, left, control with carbachol alone, right, carbachol plus 
FSFQ peptide. 

DETAILED DESCRIPTION OF THE PREFERR ED EMBODIMENTS 
The present invention discloses methods to create mode-matched peptides that 
have a high probability of binding to and/or modulating the function of target peptides, 
polypeptides, or proteins. The peptides are constructed from peptide templates derived 
from physicochemical transformations of the amino acid sequences of the target peptide, 
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polypeptide, or protein. In particular, the templates are derived from at least one of the 
following: 1) eigenvectors of the autocovariance matrices of the physicochemically 
transformed amino acid sequence of the target protein; 2) wavelet subsequence templates 
derived from a variety of wavelet transformations of the physicochemically transformed 
5 amino acid sequence of the target protein; and 3) redundant subsequence templates 
computed from the physicochemically transformed amino acid sequence of the target 
protein. 

In the peptide design methods described herein, we make new use of three 
techniques to characterize the dominant statistical wavelengths of a target polypeptide's 

1 0 physicochemical property mode (or modes) in order to generate templates for the 

construction of mode-matched peptides having a high probability of binding to and/or 
otherwise modulating, inhibiting or activating activity of the target protein, polypeptide 
or peptide (Fig. 1). The techniques are: (1) eigenfunction(s) construction from the 
convolution of the eigenvector(s) with an original data series, in which the eigenvectors) 

1 5 is determined from the autocovariance matrices of a sequentially lagged physicochemical 
property data series of the peptides, polypeptides and proteins; (2) all poles, maximum 
entropy power spectral transformation (in contrast to standard Fourier transformed power 
spectra) of the eigenfunction(s), which identifies the mode content of the 
physicochemical property data series or their eigenfunctions; and (3) discrete and 

20 continuous wavelet transformations, one-dimensional wavelet packets and multiple 
convolved wavelets (using a range of potential mother wavelets as listed above) which 
confirm the dominant statistical wavelengths of the eigenfunctions and locate them as 
phase amplitudes or absolute valued moduli in the constituent sequences. Additionally, 
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as described in detail below, the results of the wavelet transformations may locate one or 
more subsequences of the polypeptide that can serve as a wavelet subsequence template 
or an amino acid distribution source in the design of peptide or peptide-like molecules, or 
both. Similarly, a symbolic or literal template can be created directly from the 
5 subsequences so selected, or through the decomposition of single or multiple 

concatenated subsequences, to create a non-overlapping redundant subsequence template. 
The spectral modes of the polypeptides or proteins that emerge from the power spectra 
and wavelet transformations dictate the choice(s) of the eigenvector(s), alone or summed, 
which can then be used as templates for the construction of mode-matching peptides. 

1 0 The template then may be used in the manner described below to generate peptide ligands 
having a high probability of binding to and/or otherwise modulating the activity of the 
target polypeptide or protein. 

The array of potential target peptides, polypeptides, or proteins may include, 
without limitation, cell membrane receptors, nuclear membrane receptors, circulating 

1 5 receptors, enzymes, membrane and circulating transporters, membrane proteins involved 
in the translocation of viral and other infective agents into the cell, chaperonins and 
chaperonin-like proteins, monoclonal antibodies and antibody derivatives, such as Fc, 
Fab', F(ab% Fv or scFv fragments; and generally any protein, polypeptide or peptide 
involved in peptide-protein and/or protein-protein interactions. 

20 Generally, the first method of the present invention involves the linear 

decomposition of M-lagged, autocovariance matrices, Cm, constructed from the 
sequentially lagged data matrix of w . j/s of iV-length membrane proteins. M is often 
(but not always) chosen to optimize the least squares fit of the protein's leading 
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eigenfunction with its hydropathy plot, because, particularly in the case of the seven- 
transmembrane receptors, the graphs of the leading eigenfunctions closely resemble those 
of the target protein's smoothed hydropathy sequence, created by the repeated application 
of nearest neighbor averaging of the H it /=/...#. From the set of ordered eigenvalues, 
5 { v,} ,=/ ,.m of the C M , the corresponding set of ordered eigenvectors, X i( ,=/..m are computed 
and serially convolved with Hi, to form an ordered set of hydrophobic free energy 
eigenfunctions, ,=/..^, each of length N-M+l. An alternative eigenfunction 
computation, described below, results in eigenfunctions of length N. The ordered 
eigenvalue spectra generally decay quickly after the first few leading ordered values, 

10 such that most if not all of the transmembrane and peptide binding/modulating mode 
information is captured in the first few eigenvalues, i.e., {v,} w ...* although 8 < M< 25 
may be employed as required for adequate separation and resolution. 

Next, the hydrophobic mode content of the ^F/'s containing the peptide- 
binding/modulating inverse spatial frequency mode (expressed as a wavenumber, o" 1 , in 

1 5 units ofamino acids) is identified using all poles, maximum entropy power spectral 
transformations S(co) and/or wavelet transformations W(a,b). These methods revealed 
sets of statistical wavenumber matches between peptide ligands and their corresponding 
membrane receptor proteins, ranging from (o l « 2-14 amino acids across examples. 
Estimating the dominant wavenumber content of secondary eigenfunctions, ¥2, using all 

20 poles, maximum entropy power spectral transformations, and/or discrete and 

continuous wavelet and one dimensional wavelet packet transformations W(a,b), led to 
clearly resolved mode matches between peptide ligands and their receptors, and predicted 
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kinetic interactions between AGh P sequential mode-matched peptide ligands and the 
receptors. Matches as statistical patterns in AGhp modes were found between peptide 
ligands and their membrane receptors, including kappa, mu, delta and orphan opiate 
receptors, corticotropin releasing factor receptor, cholecystokinin receptor, neuropeptide 
5 Y receptor, somatostatin receptor, bombesin receptor, and neurotensin receptor. AGh P 
mode matches, such as those found between the dopamine co-localized neuropeptide 
neurotensin and the D2 dopamine membrane receptor, D2DA, and those found between 
the gastrointestinal and brain peptide cholecystokinin and the dopamine membrane 
transporter, DAT, predicted the differential binding of the pharmacologically active 

10 ligands to their respective responsive dopamine membrane receptors and, 

correspondingly, their lack of binding to the opposing, pharmacologically unresponsive 
dopamine membrane receptors. 

While the present invention is described below by employing the hydrophobic 
free energies (AGh P ) of the twenty naturally-occurring amino acids, in generating 

1 5 potential receptor binding and/or modulating peptides other quantifiable 
physicochemical properties that can order the amino acids along a particular 
physicochemical dimension of varying continuity may be used in place of the 
hydrophobic free energies. Other amino acid physicochemical properties that may be 
considered in choosing the appropriate physicochemical property include, without 

20 limitation, relative vapor pressure, relative free energy of amino acid transfer into bulk 
phases, amino acid partition coefficients in other solvent systems, diffusivity, frictional 
coefficient, aqueous cavity surface area, aqueous molar volume, partial specific volume, 
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accessible surface area, charge, mass (in daltons), volume, pK a of ionizing side chain, 
chromatographic mobility, electrophoretic mobility, chemical categorical membership 
(nonpolar aliphatic, nonpolar aromatic, polar uncharged, polar charged, basic-positively 
charged, acidic-negatively charged, sulfur-containing), structure breakers (proline, 

5 glycine), and relative occurrence in specific or groups of proteins (as percents). Other 
published properties are known to those in the art and available, for example, on the 
World Wide Web site http://www.expasy.ch. It is generally known from 
physicochemical studies that there are relatively high correlations (r = 0.6-0.8) among the 
values for the twenty naturally-occurring amino acids of free energy of transfer from 

10 aqueous to hydrophobic solvents (i.e., hydrophobic free energy), relative vapor pressure, 
aqueous cavity surface area, aqueous molar volume, partial specific volume, solvent 
accessible surface area, and other physicochemical properties. As a result, the results 
obtained from any of these quantifiable physicochemical properties would be expected to 
apply equally to the remainder of the quantifiable physicochemical properties. 

1 5 The eigenfunctions used in the eigenvector-based method are related to the 

Karhunen-Loeve, principal components and factor analysis transformations, and are 
uniquely defined in terms of an eigenvalue decomposition of each hydrophobic free 
energy data set, resulting in a set of hydrophobic free energy eigenvector-weighted 
eigenfunctions. Where available, the set of characteristic hydrophobic free energy 

20 wavelengths are isolated in the extracellular domains of transmembrane receptors. For 
example, the leading eigenfunction, V u associated with the largest eigenvalue of the 
covariance matrix of a seven-transmembrane receptor sequence locates the same 
transmembrane segments as are seen in conventional n-block averaged hydropathy plots. 
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However, unlike the case with n-block averaged hydropathy plots, the eigenfunctions 
generated by the methods of the present invention leave the remaining secondary 
hydrophobic mode (or modes) unsmoothed and available for further analyses as 
secondary eigenfunctions (i.e., ¥2, ^3,. . The eigenvectors associated with these 
5 secondary eigenfunctions may then be used as templates for the construction of mode- 
matched peptides or peptide-like molecules that then may be tested for their ability to 
bind to and modulate, activate and/or inhibit the function of the seven-transmembrane 
receptors. For other, non-seven-transmembrane receptor sequences, such as, for 
example, the human NGF receptor, the eigenvectors associated with the leading 

10 eigenfunctions may be suitable for use as peptide construction templates, since these 

hydrophobic modes are not likely to be dominated by transmembrane segments, as in the 
case of seven-transmembrane receptors. 

Alternatively, templates may be created using other methods which incorporate 
the results of discrete or continuous wavelet transformations, one-dimensional wavelet 

1 5 packet transformations or the convolution of the coefficients of two or more wavelet 
transformations. These transformations locate one or more subsequences of the target 
polypeptide that can serve as a symbolic or literal wavelet template, derived directly from 
the subsequences so selected or through the decomposition of single or multiple 
concatenated subsequences to create an eigenvector template. Still other templates may 

20 be created through the identification of symbolic or literal amino acid redundant 
subsequences in the polypeptide and peptides or peptide-like molecules known or 
believed to bind to and/or otherwise modulate the target polypeptide. 
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The methods of the present invention are described in detail below, using the 
example of hydrophobic free energy as the physicochemical property. As noted above, 
the correlations among the various physicochemical parameters enable general use of the 
methods of the present invention with other physicochemical properties, and one of 
5 ordinary skill in the art would appreciate that no undue experimentation would be 
required to perform the methods of the present invention using other physicochemical 
properties. 

A hydrophobic free energy series, H u is established for the twenty naturally- 
occurring amino acids. The values are normalized such that the reference amino acid, 
10 glycine, without a secondary structure-forming side chain, is set equal to 0.00. The 

values for Hi of each of the twenty naturally occurring amino acids cluster naturally into 
four groups, as shown in Table 1. 



Table 1 





Group 


I 


Group II 


Group 


III 


Group IV 




Amino 
Acid/Symbol 


Hi 

(kcal/mol) 


Amino 
Acid/Symbol 


Hi 

(kcal/mol) 


Amino 
Acid/Symbol 


Hi 

(kcal/mol) 


Amino 
Acid/Symbol 


H t 

(kcal/mol) 




tryptophan/W 


3.77 


cysteine/C 


1.52 


alanine/A 


0.87 


serine/S 


0.07 




tyrosine/Y 


2.76 


methionine/M 


1.67 


aspartate/D 


0.66 


threonine/T 


0.07 




phenyalanine/F 


2.87 


valine/V 


1.87 


histidine/H 


0.87 


glycine/G 


0.00 




isoleucine/I 


3.15 


lysine/K 


1.64 


arginine/R 


0.85 


glutamine/Q 


0.00 




proline/P 


2.77 


leucine/L 


2.17 


glutamate/E 


0.67 


asparagine/N 


0.09 



1 5 The set of hydrophobic free energy values naturally clusters into four 

discontinuous groups, with two exceptions. Proline (P), though having a value of 2.77 
kcal/mol which places it in the highest hydrophobicity group, acts as a secondary 
structure breaker, due to its rigid constraints on rotation about the N — Ccc bond and 
absence of an amide hydrogen for resonance stabilization of its peptide bond or 
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participation in carbonyl-imino H-bonding. Consequently, proline has unusual hydrogen 
binding inclinations and <€ breaks" the continuity of one-dimensional hydrophobic waves 
in the same way as its nucleotide complement partner in the lowest hydrophobicity group, 
glycine. Therefore, proline is assigned to the lowest hydrophobicity group with glycine 
5 and is given the same value (see Table 2). In addition, leucine has many of the properties 
of the highest hydrophobicity group and is assigned to that group in place of proline. 
Therefore, the twenty naturally occurring amino acids are divided on the basis of the 
hydrophobic free energy values into four hydrophobicity groups consisting of the 
following amino acids: Group I (highest hydrophobicity): L,W,Y,F,I; Group II (second 
10 highest hydrophobicity): C,M,V,K; Group III (third highest (second lowest) 

hydrophobicity): A,D,H,R,E; and Group IV (lowest hydrophobicity): S,T,G,Q,N,P. 
These groupings are shown in Table 2. 



Table 2 



Group 


I 


Group II 


Group ffl 


Group IV 


Amino 
Acid/Symbol 


H t 

(kcal/mol) 


Amino 
Acid/Symbol 


Hi 

(kcal/mol) 


Amino 
Acid/Symbol 


Hi 

(kcal/mol) 


Amino 
Acid/Symbol 


Ht 

(kcal/mol) 


leucine/L 


2.17 


cysteine/C 


1.52 


alanine/A 


0.87 


serine/S 


0.07 


tryptophan/W 


3.77 


methionine^ 


1.67 


aspartate/D 


0.66 


threonine/T 


0.07 


tyrosine/Y 


2.76 


valine/V 


1.87 


histidine/H 


0.87 


glycine/G 


0.00 


phenyalanine/F 


2.87 


lysine/K 


1.64 


arginine/R 


0.85 


glutamine/Q 


0.00 


isoleucine/I 


3.15 






glutamate/E 


0.67 


asparagine/N 


0.09 














proline/P 


0.00 



1 5 The natural division of H into four sets of four to six amino acids each (Tables 1 

and 2) is used in assignment of amino acids to the four-partitioned eigenvector templates 
used in the construction of new candidate peptide ligands, while the values of H t in Table 
2 are used in the transformation of the amino acid sequence of the receptor into a real 
number AGh P series, as described below. It will be apparent to one of skill in the art that 
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other groupings are potentially appropriate and that as other physicochemical properties 
are employed, the amino acids may group differently. 

Each target polypeptide having an amino acid sequence of length N 9 comprised of 
amino acids Aj,A 2y .A N may be represented as a sequence of hydrophobic free energy 
5 values H h H 2 , . . ~H N , where Hi represents the hydrophobic free energy value of amino 
acid 4- in the /-th place in the amino acid sequence, using the H t values listed in Table 2 
above. Each target polypeptide sequence, H h H 2 , . . .H N , is transformed first into a 
sequentially lagged data matrix, then into an autocovariance matrix, and finally 
decomposed into a set of orthogonal functions. 

10 From the data column vectors (T= transpose) F/= {H } , H 2 ..„ H n . M \ V 2 (H 2 , 

H 3 ,...,H*mh), - V M T = (H m ,Hmh H n ) and where K = the sequence averaged 

dyadic product, HiHi is used to obtain the autocovariance matrix, a Hermitean MxM 
matrix, C M = 1/K { HiHi } . M is sometimes chosen to minimize the least squares error of 
the protein's leading eigenfunction,^!, with their hydropathy plots resulting from the 

1 5 standard technique of nearest-neighbor averaging. As such, values for M are often in the 
range of about 10 to about 20. 

The eigenvalues, { V /} M and the associated eigenvectors, Xi(j), of Cm, are 
calculated wherein / = 1 . . M and labels the eigenvector, andy = 1 . . M and refers to the 
yth component of the eigenvector^^. The eigenvalues {v*}^ are ordered from largest 

20 to smallest, as are the corresponding eigenvectors Xi(j). The ordered Xi(j) are then used 
as multiplicative "weights" to transform the H lt H 2 , ... f H N into M statistically weighted 
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eigenfiinctions, ^(j), where i = 1 . . M labels the eigenfunction and 7=1 . . .N-M indexes its 
yth component. The ^{j), for j-k+\ > 0, are given by 

M 

Alternatively, N length ^(j), for j > 0, are given by 

Here Hi is the first hydrophobic free energy value in the sequence. Intuitively, 
Cm scans for hydrophobic modes across a range of autocorrelation lengths from 1 to M, 
the range of the lags in the autoco variance matrices. Because Cjtf is real, symmetric (Hy 
=H 9 ) and normal (CmC m T =C m T C m \ its { m} (W ... m are real, non-negative and distinct, and 
its associated eigenvectors, X$) 9 constitute a natural basis for orthonormal projections on 
H h H 2 ,...JIn- The set of can be regarded as orthonormally decomposed sequences 
of eigenvector-weighted, moving average values. 

The eigenfiinctions may be shortened with respect to the receptor sequences by 
the number of lags Mused to construct the covariance matrix. The leading eigenfunction 
representing the transmembrane segments of receptor proteins is designated as *¥ T 9 the 
secondary eigenfiinctions containing the peptide-binding/modulating receptor mode or 
modes as *F*, and the leading peptide or peptide-like molecule ligand eigenfunction as 
(when the peptide or peptide-like molecule is long enough to permit its construction). 
The eigenvectors serve as weights to generate orthonormally decomposed sequences of 
moving average values with the potential for finer resolution of mode or modes than that 
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possible in the moving average graph of the hydropathy plot or Fourier transformation of 
the undecomposed data series. 

In the computation of maximum entropy power spectral transformations, S(cd), the 
a* coefficients are calculated directly from the H t or series, and represent the average 
5 over Hi separated by k residues or values in the relevant sequence such that = 

{H(i)H(i + k)) = — ?— k H ®H(i + *) for N ' M + l P oints in the case of the ^ 

ftj k 

Where z = e i6) , the conventional Fourier power spectral transformation is inverted such 

1 



that poles replace the zeros of the usual expansion; i.e., in S((o) = 



N-M+l 

1+ Xajtz 

k=\ 



, where 



the denominator is a minimum, S(®) will have peaks. It can be shown, using the method 
1 0 of Lagrange multipliers, that extending beyond the known a* s for k = -M. . M into a 
Gaussian process maximizes the entropy, H 9 ofS((o\ H = j\nS(co)d6) in the all poles 

power spectral transformation. Here, k is the number of poles chosen for examination 
and is usually (but not always) held to < 8 for receptor eigenfunctions derived from 
receptors having sequence lengths of several hundred amino acids to avoid "splitting" S 
1 5 (co) into spurious modes. In the all poles maximum entropy power spectral 

transformation, 5(g)) is much like an autoregressive, maximum-likelihood spectral 
estimate in that it is not mode-dependent, but is derived directly from the data of H u = L „n 
and *F, and behaves like a filter that may yield the one or two leading poles of discrete 
hydrophobic variation in the hydrophobic free energy eigenfunction. 

20 
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Wavelet Transformations of Hydrophobic Free Energy Functions 

Whereas the of the protein's leading and its ligand's *¥ L locate the 

conjectured binding/modulating mode or modes in AGh P wavelength space, their 
sequence position is lost. In contrast, wavelet transformations yield sequence and 
wavelength information simultaneously. Discrete wavelet techniques allow cutting 
smooth windows of differing lengths while preserving orthogonality during pattern 
identification in W. 

Haar, Trigonometric, Meyer, Daubechies, Gabor, Battle-Lemarie, Biorthogonal, 
Coifinan, Grossmann, Morlet, Mexican Hat and other mother wavelet families may be 
used in wavelet transformations that depict specific proteins as a signatory sequence of 
hierarchical modules. These functions are called wavelets because they have a local 
oscillatory form, so that, unlike the sinusoidal waves of Fourier transformation, they 
decay as H -> ». There are a wide variety of choices of "mother wavelets" which are 
systematically dilated, translated and then composed with the original sequence. 

With respect to the hierarchical scaling characteristics, unlike the Fourier 
transform which sacrifices location for knowledge of characteristic wave numbers, the 
wavelet transformation is well suited to study regions of non-random autocorrelation 
which appear intermittently across a sequence and with hierarchies of scale. This is 
exemplified in proteins by the typical patterns of alternating helices, strands and loops as 
localized coherent structures along longer wavelengths of intermittent patterns of larger 
autocorrelated sequential structures, such as helical barrels and sheets. These, in turn, ar 
components of still larger autocorrelated sequences in the form of protein domains. 
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With respect to hydrophobic free energy sequences, we have found that the 
Dubechies wavelets, and in particular its simplest member, the Haar wavelet, are usually 
better suited for locating structures in sequence space, while the Morlet, Meyer and 
Mexican Hat wavelets are best for indexing sequential structures in dilate space. The 
approach using the Morlet mother wavelet is presented here. However, it will be 
understood by those of skill in the art that other mother wavelets could also be employed, 
as desired. The wavelet method of locating, describing mode relevant subsequence and 
constructing wavelet subsequence templates from which to design peptides for binding, 
modulation, activation and/or inhibition of a target polypeptide/protein is has not been 
previously described and is unique to the present invention. 

Assuming the protein structural organization that was first suggested by 
Linderstrem-Lang and assuming, for example, 64 dilate divisions are in the wavelet 
graph, some or all of the following kinds of information are available from the Morlet 
wavelet transformations of an undecomposed H t . First, at relatively small scales, the 
sequence locations and fundamental sequential hydrophobic inverse spatial frequencies 
or wavenumbers of the protein's characteristic secondary structures can be determined. 
For example, a-helices contain from 3.2 to 3.7 amino acids per hydrophobic free energy 
rotation (« 24-30 dd), while p-strands have rotation numbers which may range between 
2.2 to 2.6 amino acids (« 5 to 15 dd). Second, at intermediate scales, the characteristic 
sequence sizes and locations of singular, hierarchical, secondary structures can be 
assessed. For example, although there is considerable variability, individual helices in 
helical bundles generally average in the range of 7 to 15 residues in length (« 48 to 55 dd) 
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and p-strands in sheets or barrels may range from 4 to 8 residues 32 to 45 dd). Third, 
at the next largest scale, the multiresolution capacity of W{a f b) may be exploited to locate 
another kind of sequence similarity characteristic of the multiscale, hydrophobic 
sequence content of the longer and shorter loops (called "random coils"), which serve as 
5 transitions between more dilate localized secondary modules of helices or sheets. These 
random coils range generally from 2 to 16 residues, although they can be longer. Lastly, 
the modular maxima at the largest scales (> 60 dd) are relatively long hierarchical 
hydrophobic domains of 40 to 50 amino acids, or more. 

The complex Morlet continuous wavelet transformation, W{a,b\ of a protein's 
10 undecomposed Hi is obtained by dilating (i/a) and translating (i/b) the analyzing wavelet, 
w. With b representing distance translated down the sequence and a the "scales" or 
"dilates" as sequential radian frequencies or wavenumbers of w, the "mother wavelet", 

wavelet transformations, W(a f b) = (1/ 4a) [H(i)w(- — -)di may be useful in conserving 

o a 

both wavelengths and locations for structural prediction using Hi in polypeptides and 

1 5 proteins. For w we chose a member of the family of continuous, symmetric, » zero mean, 

infinitely regular and differentiable, modulated Gaussian Morlet wavelets - 

1 - x 2 

w(x) - — exp( ) cxp(2mfic) . 

2n 2 

Even though this and most of its other applications involve real numbered series, 
the Morlet continuous wavelet transformation W{a,b) is complex. As such, it has real 
20 (modulus) and imaginary (phase) parts. In categorizing proteins into structural families, 
the physicochemical features (i.e., hydrophobic free energies or other amino acid physical 
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properties listed above) of the sequence locations, wavenumbers and hierarchically 
scaling transitions are of interest. Both the phase and modulus plots are suited to the 
detection and location of such features. 

Intuitively, the usual three-dimensional wavelet space (not shown) exploits 64 
5 dilate divisions, dd, related to mother wavelengths, w, as a nonlinear function, w =J{dd) 

1 b To prevent aliasing, the shortest m = 1/0.5 = 2 amino acids, which is 

0.5-(dtfX^j) 
64 

graphed at the bottom end of the y-axis, with/^^l/O = oo at the top end. The position 
on the x-axis indexes sequence location; the y-axis indicates the relative dilation of w(x) 
(composed with Hi) in dilate divisions. The modular amplitudes of the wavelet 

10 transformations may be graphed as gray-scale shaded, with relative maxima being lighter 
and relative minima being darker in shading. These absolute amplitudes within each of 
the 64 dilate ranges were normalized to 100 % ("coloration by scale'*). This choice of 
"by scale" versus "across scale" color coding of modular amplitudes does not portray the 
relative dominance of structures across all dilate ranges (which results in the loss of 

15 wavelet structural detail), but rather outlines the relative amplitudes of modular patterns 
and their locations at each dilate range. A variety of graphing techniques including color 
coding, gray scale, contour and other ways of indicating moduli and/or amplitudes may 
be employed, as determined by the particular global polypeptide property that is being 
addressed. 

20 The wavelet transformation method transforms a one-dimensional Hi series into a 

two-dimensional wavelet space, resulting in informational redundancy that is inherent in 
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the wavelet transformation technique. Potentially artifactual autocorrelations due to the 
redundancies can be defined in terms of their average over the entire sequence of 
observables. It is known, for example, that continuous wavelet graphs of random series 
can manifest patches of correlated regions which decrease with increasing scale and have 
5 their origins in the wavelet of the transform itself. In light of this problem, the Morlet or 
other wavelength graphs of the eigenfunctions, as opposed to those of the undecomposed 
sequences, may be used to seek additional information in support of the origins in the 
data of the structural features of the wavelet graphs. 

Wavelet transformations of the receptor and ligand eigenfunctions generate 
10 wavelet graphs, W* and FF 6 . Wavelet transformation, W(a,b) of the receptor 

eigenfunction W K is accomplished by decomposing the eigenfunction values into 
translated W(n)-+ W{n~b) and scaled W(n)-> W(n/a) versions of the mother wavelet, w, a 

waveform having an average value of 0 ^J* w(n)dn = 0 j , of finite length, arbitrary 

' i — b 

regularity and symmetry, and which is composed as W(a,b) = (1 / 4a) \H{i)w{ )di , 

0 ^ 

1 5 as above for the undecomposed if/. Similarly, wavelet transformation of the ligand 
eigenfunction ^ is accomplished in the same manner, by decomposing the T L values 
into translated W{n)-± W{n-b) and scaled W(n)~* W(n/a) versions of the mother wavelet, 
w. 

Because wavelet transforms preserve sequence position information of the 
20 statistical modes 1 occurrences, the results of any of the variety of wavelet transformations 
locate one or more subsequences of the polypeptide that can serve as amino acid 



57 



PATENT 

EXPRESS MAIL NO. EL557885546US 
distribution sources in the design of peptide or peptide-like molecules. The distribution 
of amino acids within these subsequences can be employed as a guide in the selection of 
particular amino acids within the physicochemical group of the peptide template, as 
further discussed below. Similarly, a symbolic or literal template can be created directly 
5 from the amino acid subsequences corresponding to the physicochemical subsequence or 
subsequences so selected or through the decomposition of single or multiple concatenated 
subsequences to create an eigenvector template. 

While the peptides or peptide-like molecules produced by this method almost 
always share the maximum entropy power spectral modes of their eigenvector template, 

10 it is sometimes the case, particularly when the eigenvector template is multimodal, that a 
mode evident in the maximum entropy power spectrum and wavelet transformations of 
the eigenfunction or eigenfunctions of interest is not evident in the maximum entropy 
power spectral transformation of the associated eigenvector or eigenvectors, their 
template or the peptides produced from the template. Often the spectrally invisible mode 

15 has the longer wavelength of multiple modes, and when this is the case, the mode is often 
detectable as an amplitude-modulated wave in the eigenvector, its template or the 
peptides produced from the template. This may result from the short length of the 
eigenvector, its template and the peptides produced from the template and the statistical 
nature of the power spectral transformation. The eigenvector, its template and the 

20 peptides produced from the template are still considered to be mode-matched to the 
polypeptide, as they contain physicochemical amplitude variations on the mode of 
interest. 
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Wavelet Packet Transformations 

Wavelet packet analysis may also be used in the identification, localization and 
characterization of physicochemical modes and mode relevant subsequences and the 
creation of wavelet subsequence templates. Wavelet packet analysis uses the same set of 

5 mother wavelets listed above, but generalizes the technique, allowing a range of 
representations of the decomposed sequence. In one-dimensional wavelet packet 
analysis, the physicochemical series, S, is decomposed into its gross and fine scale 
variation, then each of the resulting gross scale, G, (approximation) and fine scale, F, 
(detail) series are again decomposed into gross and fine scales. This process is repeated 

1 0 an arbitrary number of times, p, resulting in a binary tree of sequences with p levels as 
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20 



for p=2. The original physicochemical series can then be represented as an expansion of 
the wavelet packet atoms, each of which is a waveform, e.g., S=Gi+GF 2 +FF 2 . As p 
increases and trees get more complex, the number of such possible representations is 
obviously large. To select among these representations of the physicochemical series we 
tploy one of two entropy threshold criteria: Shannon (i.e., # 2 log(# 2 ) ) and Stein's 



em- 



Unbiased Risk Estimate (SURE) (i.e., V21og e (nlog 2 («)) , where n equals the number of 
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points in the physicochemical series). With these criteria we produce "best level" and 
"best tree" representations, with which we can compare the physicochemical attributes of 
two or more physicochemical series. 

Wavelet packets are relatively easy to compute when using orthogonal mother 
5 wavelets. Starting with two filters of length N corresponding to the wavelet, h(n) and 
g(n), the reversed version of the low-pass decomposition filter and the high-pass 

decomposition filter are divided by 4l respectively. Then we define the system of 
functions W n (x\ (n=0,l,2...) as, 

2N-1 

10 and 

2JV-1 

W 2n+l {x)^g(k)W n {2x-k) 

k=0 

where W 0 (x) is the scaling function and W\(x) is the wavelet function. 

Starting from the functions W n {x), neN, we consider the family of analyzing 
functions W M>k (x)=2~ J/2 W n (2~ J x-k) , where n e AT and jjc are nonnegative integers.; can 

15 be considered a scale parameter and k can be interpreted as the sequence localization 
parameter. W n (x) oscillates approximately n times. For fixed j and £, Wj^k assesses 
fluctuations of the physicochemical sequence around the position 2 J • k at the scale 2~ J 
across frequencies/wavenumbers for the accessible values of n. For some basis functions, 
the naturally ^-ordered functions must be reordered so that the number of zero crossings 

20 of the wavelet increases monotonically with the order of the function. 
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The set of functions W/^x) is the (j,n) wavelet packet, which when j,n are 
positive integers and k has an integer value, are organized in tree structures. Each node 
of the tree is of the form 

A 

Because {(Wj + j f2 „), (Wj+]j n +i)} is an orthogonal basis of the space spanned by 
W Jt n 9 the leaves of every connected binary subtree of the wavelet packet tree correspond 
to an orthogonal basis of the initial space. For our physicochemical sequences, each 
wavelet packet basis will provide an exact reconstruction but with a specific spatial 
frequency subband coding. As a result, a physicochemical series of length N= 2 L can be 
expanded in at most 2 N ways with a binary tree of depth L. 

As these can be unmanageably large numbers, we choose optimal representations 
through the application of the two entropy criteria listed above, i.e., Shannon and Stein's 
Unbiased Risk Estimate, although other criteria could be employed. Other entropy-based 
criteria usable in the wavelet packet transformations can include the logarithm of the 
"energies" entropy (i.e., £ log(Hf), with the convention that log(0)=0), topological 

entropy estimate for a finite series (Le., the asymptotic growth rate of the trace of the 
recursively exponentiated transfer matrix of each subband), and a fixed entropy 
threshold. Because they are well suited to quantifying additivity type properties, produce 
efficient searches in binary tree structures, and describe information carrying properties 
of the subbands, we favor entropy-based criteria. 
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In each case, we compute the entropy of the original physicochemical series, then 
we split the series using the chosen wavelet and recompute the entropy of each resulting 
piece. If the sum of the entropies of the pieces at a given level is less than the sum of the 
entropies of the preceding level, the split is considered to be informative. By this method, 
applied exhaustively to all possible additive representations, entropy-minimizing best 
level and best tree representations can be defined. These graphs are frequency-ordered 
(i.e., subband graphs are arranged from those representing lowest to those representing 
highest frequencies) so as to be maximally interpretable. A variety of graphing 
techniques, including "by scale" and "across scale" color coding, gray scale, contour and 
other ways of indicating coefficient values may be employed. 

Intersection of Two or More Wavelet Coefficient Arrays 

The intersection of two or more wavelet coefficient arrays may also be used in the 
identification, localization and characterization of physicochemical modes and mode 
relevant subsequences and the creation of wavelet subsequence templates. Various 
wavelet techniques are differentially suited to the assessment of specific aspects of the 
physicochemical protein, polypeptide and peptide or peptide-like molecule series. For 
example, as noted above, in discrete or continuous wavelet analysis, Haar mother 
wavelets are particularly suited to localizing coefficients in sequence space, while Meyer, 
Morlet and Mexican Hat mother wavelets are better suited to dilate space localization. 
To derive more information in a single representation, and if the matrices of coefficients 
are of the same order and derived from analyses of the same physicochemical series, we 
generally apply highpass filters to each wavelet coefficient matrix and then compute their 
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cell-wise intersection. A nonzero cell, Ag, in each and all constituent matrices results in 
a nonzero corresponding cell in the intersection matrix, By, that takes a value equal to the 
average or median of the values of the corresponding cells in the constituent matrices. 
Constituent wavelet coefficient arrays can result from the use of discrete or continuous 
5 wavelet transforms or wavelet packet analysis, and from any of the above listed mother 
wavelets, provided the above conditions are met. The intersection matrix serves to 
evaluate the wavelength and dominant position or positions of physicochemical modes, 
and also as a method by which to identify one or more amino acid subsequences in the 
analyzed polypeptide or peptide that are associated with mode-relevant binding and/or 

1 0 modulation. The subsequence or subsequences so identified may be employed 
individually or together as a source for amino acid probabilities in the creation of 
peptides or peptide-like molecules. The amino acid or corresponding physicochemical 
subsequence or subsequences may be used directly or in a coded form as a template for 
the design of peptides or peptide-like molecules that will bind the polypeptide or peptide 

1 5 on which the analysis was based. 

Construction of Peptides by Assignment of 
Amino Acids to An Eigenvector Template 

20 The sequential eigenstructures of the transformations described above may be 

used to design de novo new peptides that may bind to and/or otherwise modulate and 

have an influence on various protein or polypeptide activities. To construct new peptide 

ligands, the sequential H t (or other physicochemical properties, as above) values of the 

receptor are normalized and partitioned. Amino acid assignment is dictated by the mode- 
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relevant eigenvector or eigenvector-based template, and is consistent with membership in 
one of the natural divisions dictated by the physicochemical property, e.g. the four natural 
divisions of the naturally-occurring amino acid's AGh P values. Furthermore, amino acid 
assignment may be weighted by any desired means known to those in the art, such as by 
5 the amino acid distribution found in a particular amino acid pool or by accounting for 
known effects of directed mutations or segment replacements. 

Peptide construction from the distinct spectral signature eigenvector-based 
template begins with the selection of the appropriate eigenvector (or eigenvectors), based 
on their eigenvalues and the maximum entropy power spectral mode or modes of the 

1 0 associated eigenfunction or eigenfunctions to be represented in the eigenvector template, 
X temp . The y-axis of the graph of X temp is divided into a number of segments, 
corresponding to the range of AGh P values of each of the various groups of the twenty 
essential amino acids listed above in Table 1 or Table 2. The index of the eigenvector 
(graphed on the x-axis ofX temp ) maybe any value between 1 and M, and is chosen based 

1 5 on the relevant eigenfunctions that the all poles power spectrum and/or the wavelet 
transformation have shown contain the receptor's ligand-matching signatory mode or 
modes. For example, in the cases of the seven-transmembrane receptor superfamily 
members, the first eigenfunction (i = 1) resembles the moving average hydropathy plot, 
and it is the second (and sometimes additionally a higher eigenfunction) that provides the 

20 distinct spectral signature of the protein that may act as the template for the construction 
of the mode-matched peptide. In the cases of the single transmembrane tyrosine kinase- 
coupled receptors, and other receptors with a single transmembrane sequence (and other 
protein families listed above), as well as other proteins, such as transporters, enzymes and 
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chaperones, the first eigenfunction (and again, sometimes additionally a higher 
eigenfunction) may contain useful spectral signatures. The ordered eigenvalue spectra 
generally decay quickly after the first few leading ordered values, such that most if not all 
of the transmembrane and peptide binding/modulating mode or modes information is 
5 captured in the first few eigenvalues, i.e., {v,} ... 4 , though 8 < M < 25 may be employed 
for adequate separation and resolution. 

With respect to the substitution process in the M-length eigenvector template Xtemp 
associated with the eigenfunction or eigenfunctions of interest, the sequence of values in 
the ^(vector position)-><vector position) of X temp are plotted, followed by partitioning of 

1 0 the occupied region of the y axis into the desired number of parts. While the hydrophobic 
values of the twenty naturally-occurring amino acids naturally partition into four equal 
parts (Table 1 and Table 2), the hydrophobic values may also be partitioned into a lesser 
or greater number of parts, and the partitions may or may not be equal. Furthermore, 
when other physicochemical properties are used, another number of partitions may be 

1 5 desirable. In the case of hydrophobic free energies, the top region of the partitioned 

eigenvector template graph is mapped to the highest hydrophobicity (i.e., Group I) amino 
acids, the next region to the second highest hydrophobicity amino acids (i.e., Group IT), 
etc. down to the lowest hydrophobicity amino acids in the lowest region. Starting at the 
first of the M points of X temp , the amino acid hydrophobicity group to which this point 

20 belongs is determined. Then, a member of the amino acids in this group (from the chosen 
amino acid pool) is randomly assigned to this point. The process then is repeated for the 
remainder of the points in the eigenvector template to generate an M-length peptide 
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which is considered mode-matched to the receptor. The process may be repeated as often 
as desired to generate a large number of eigenvector template-defined candidate peptides. 

Multiple eigenvectors derived from the same receptor (e.g., X\ and X 2 ), each with 
distinct spectral properties in the associated eigenfunctions may also be used in 
5 combination to generate candidate peptides. In such a case, it is important to preserve 
multiple aspects of the receptor's eigenfunction mode signature. Accordingly, an 
eigenvector template vector Q of length M is formed. Vector Q is the eigenvalue (v)- 
weighted sum of the eigenvectors (X) from which the eigenfunctions are derived. That is, 
Q(j) = v\X\ + v 2 X 2 . This is possible due to the linear additivity of eigenvectors and their 
10 eigenvalue weights. The candidate peptides then are generated as described above, using 
Q(j) in place of the single eigenvector in the assignment of the amino acids from the four 
amino acid groups. It will be obvious to those of skill in the art that other transformations 
and composites of multiple eigenvectors can be employed to form M-length eigenvector 
templates derived from two or more eigenvectors, as desired. 

15 

Construction of Peptides by Assignment of Amino Acids Based on Mode-Matching 
to Wavelet Identified Polypeptide Subsequences 

Like the sequential eigenstructures described above, the results of the variety of 
wavelet transformations described above may be used to design de novo peptides that 
20 may bind to and/or otherwise modulate and have an influence on various protein or 

polypeptide activities. The wavelet-derived subsequence template, S te mp> is produced by 
first performing discrete or continuous wavelet transformations, wavelet packet 
transformations or multiple convolved wavelet transformations on a polypeptide 
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physicochemical series and on the physicochemical series of a peptide or peptide-like 
molecule known or suspected to bind the polypeptide. Modes of physicochemical 
fluctuation are assessed to identify the mode or modes of interest, generally as a mode or 
modes shared by the polypeptide and the peptide under consideration. Once this mode or 
5 modes is identified, it can be localized in the sequence of the polypeptide by selecting an 
interval around wavelet coefficient peaks in the dilate subband or subbands that 
correspond to that wavelength. These sequence intervals are then used to select the 
corresponding sequences of amino acids in the primary polypeptide series. Amino acid 
subsequences are then coded into group membership on the basis of a physicochemical 

1 0 property and its grouping scheme, and this coded sequence acts as a template for the de 
novo generation of new peptides, as above. 

To construct new peptides, the sequential physicochemical values of the 
polypeptide or protein and peptides known to bind it, if such exist, -are normalized and 
partitioned. Shared physicochemical mode or modes or mode(s) of interest are identified 

15 in the wavelet graphs. The sequence interval at which a mode is dominant in the 

polypeptide is identified and this subsequence of 100 amino acids or less in length forms 
a template. Amino acid assignment is dictated by the mode-relevant subsequence-based 
template, and is consistent with membership in one of the natural divisions dictated by 
the physicochemical property, e.g. the four natural divisions of amino acid AGh P . 

20 Furthermore, amino acid assignment may be weighted by any desired means known to 
those in the art, such as by the amino acid distribution found in a particular amino acid 
pool, or by accounting for known effects of directed mutations or segment replacements, 
as described below. 
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The subsequence-based template, S te mp> is graphed so that the y-axis of the graph 
of Stemp is divided into a number of segments corresponding to the group memberships of 
the essential amino acids, listed above in Table 1 or Table 2. The chosen polypeptide 
amino acid subsequence may be any contiguous interval of 100 amino acids or less in 
5 length, and is chosen based on the colocalization in the dilate space shown contain the 
receptor's ligand-matching signatory mode or modes and the sequence space 
corresponding to the chosen interval 

With respect to the substitution process in the subsequence template S te mp 
associated with the subsequence or subsequences of interest, the sequence of values in the 

10 ^(vector position)-><physicochemical group membership) of S t€ mp are plotted. While the 
hydrophobic values of the twenty naturally-occurring amino acids naturally partition into 
four equal parts (Table 1 and Table 2), the physicochemical values associated with the 
amino acids may also be partitioned into a lesser or greater number of parts, and the 
partitions may or may not be equal. Furthermore, another number of partitions may be 

1 5 desirable. In the case of hydrophobic free energies, the top region of the partitioned 

template graph is mapped to the highest hydrophobicity (i.e., Group I) amino acids, the 
next region to the second highest hydrophobicity amino acids (i.e., Group II), etc. down 
to the lowest hydrophobicity amino acids in the lowest region. Starting at the first point 
of Stemp, the physicochemical value of this point is related to the appropriate 

20 hydrophobicity group. Then, a member of the amino acids in this group (from the chosen 
amino acid pool) is randomly assigned to this point. The process then is repeated for the 
remainder of the points in the template to generate a peptide which is considered mode- 
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matched to the receptor. The process may be repeated as often as desired to generate a 
large number of subsequence template-defined candidate peptides. 

Construction of Peptides by Assignment of Amino Acids Based on 
5 Redundant Polypeptide Amino Acid Subsequences 

An alternative template based on symbolic dynamics may also be used to design 
de novo peptides that may bind to and/or otherwise modulate and have an influence on 
various protein or polypeptide activities. A redundant subsequence template, R te mp> 
results from the evaluation of the symbolically-coded amino acid sequence of a target 

10 polypeptide and/or protein. The polypeptide amino acid sequence of length N is either 
retained as a string vector of amino acid one-letter representations or is transformed into a 
symbol sequence by replacing each amino acid with a value representing its group 
membership associated with a physicochemical property and a grouping scheme. In 
either case, the //length sequence, D if i=i,2 v ..M is treated as a string vector and examined 

1 5 for redundant substrings. 

Starting at the first points of the sequence, a search is made for the largest 
possible repeated substring, of length N/2, that is, points D [1,2,.. JV/2]. Next the search 
sequence size is reduced by 1, and starting again at the first point, the first N/2 - 1 
characters are assigned as a search string and all identical non-overlapping substrings are 

20 identified as the algorithm looks down the entire N length series. When this search is 

complete, the search string is reassigned as points corresponding to points D [2,3 . . N/2] 
and all non-overlapping substrings identical to the search string are identified as the 
algorithm looks down the entire N length series from Dm+i to A^. When this search is 
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complete, the search string is reassigned as points corresponding to points D [3,4..,, 
M2+1], and so on. When all possible non-overlapping redundant substrings of a given 
length have been identified, the search string length is reduced by one and the search is 
resumed. This recursive search terminates when the search string is only one character 
5 long. Redundant substrings of three or more characters must be repeated at least twice to 
be considered, while substrings of two characters must be repeated at least three times. 

All non-overlapping substrings (i.e., those with at least two distinct occurrences in 
Di) are saved and displayed with their corresponding frequencies of occurrence and 
starting positions in the A . Rtemp may be composed of a single or multiple redundant 

10 substrings so identified. When multiple substrings, or redundant substrings, are 

employed the multiple substrings are concatenated to form Rtemp* Preference is generally 
given to long subsequences in the creation ofR temp . However, the choice of redundant 
substring or substrings represented in the R te mp may be modified by knowledge of the 
results of studies of point mutations and/or peptide segment exchanges that affect 

15 binding/and or activity of ligands for the receptor, and/or specific subsequence 
physicochemical attributes from the literature. 

With respect to the substitution process in the redundant subsequence template 
Rtemp associated with the subsequence or subsequences of interest, the sequence of values 
in the ^(vector position)-,y(physicochemical group membership) are plotted. While the 

20 hydrophobic values of the twenty naturally-occurring amino acids naturally partition into 
four equal parts (Table 1 and Table 2), these or other physicochemical values associated 
with the amino acids may also be partitioned into a lesser or greater number of parts, and 
the partitions may or may not be equal. Furthermore, another number of partitions may 
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be desirable. In the case of hydrophobic free energies, the top region of the partitioned 
template graph is mapped to the highest hydrophobicity (i.e., Group I) amino acids, the 
next region to the second highest hydrophobicity amino acids (i.e., Group II), etc. down 
to the lowest hydrophobicity amino acids in the lowest region. Starting at the first point 
5 of Rtemp, the physicochemical value of this point is related to the appropriate 

hydrophobicity group. Then, a member of the amino acids in this group (from the chosen 
amino acid pool) is randomly assigned to this point. The process then is repeated for the 
remainder of the points in the template to generate a peptide which is considered mode or 
modes-matched to the receptor. The process may be repeated as often as desired to 
10 generate a large number of subsequence template-defined candidate peptides. 

As an example of redundant substring template generation, consider the following 
short amino acid sequence retained as a string vector of amino acid one letter 
representations; 

AIRCKSMLRYGHAMQLREWVCCMHAMQVYRLM 
15 If we chose to apply the template-generating algorithm directly to this series, the 

search algorithm would begin by looking for two copies of the first half of the series, 
AIRCKSMLRYGHAMQL. Next it would assess the starting positions and frequency of 
occurrence of the substring from which the last amino acid, L has been dropped, i.e., 
AIRCKSMLRYGHAMQ, and so on, looking at each possible substring in the first half of 
20 the sequence. The algorithm finds one redundant substring, HAMQ, occurring twice 

starting at positions 12 and 24. A generalization of this method also allows for the search 
of substrings that are both "backward" and "forward" in orientation in the original 
sequence. Such a search of our example string also turns up the twice repeated substring 
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MLRY, appearing at starting position 7 in a "forward" orientation and at starting position 
29 in a "backward" orientation. Our R temp might then equal one or both of these specific 
amino acid substrings in some order and orientation. 

Transforming the amino acid sequence into a symbolic vector in which each point 
5 represents the hydrophobic free energy group membership of the corresponding amino 
acid sequence, we get: 31322422314332423312222332421322. A search of this string 
for redundant substrings yields: 33242 (which appears twice starting at positions 12 and 
24), 1322 (which appears twice at starting positions 2 and 29, and corresponds to the 
MLRY sequence described above) and 22 (appearing three times that do not overlap the 
10 longer coded subsequences at positions 7, 20 and 22). Our R temp might then include one 
or any combination of these substrings representing hydrophobic free energy groups. 
Examples of appropriate sample subsequence templates might then include 
33242221322, 13222233242, and 2233242221322, among others. 



1 5 Reduction of the Number of Potential Mode-Matched Peptide Candidates 

The large number of potential candidate peptides generated in this fashion can be 
reduced in a number of ways. First, all poles power spectral analyses or wavelet 
transformations of the peptides may be performed to determine those peptides having the 
best mode-match to the receptor. In addition, the probability of occurrence of the amino 
20 acid members of each of the four AG hp groups in the general amino acid pools available 
to the particular organ or organism may be determined, and the assignment of the amino 
acids may be weighted accordingly. Finally, the results of studies of point mutations 
and/or peptide segment exchanges that affect binding and/or activity of ligands for the 
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receptor from the literature may lead to empirical attempts to optimize the sequences of 
the candidate peptides. The distributions of amino acids from which random selection by 
partition memberships may be made include the amino acid compositions of relevant 
proteins, free amino acid pools from brain, liver and/or other organs, bound and/or free 
5 amino acids in plasma and/or spinal fluid, extracellular, intracellular and/or other free 
amino acid pools, or may be derived from a subsequence or subsequences of amino acids 
located through the application of wavelet transformations or through the calculation of 
redundant substrings of amino acids. Furthermore, any combination of these procedures 
may be employed to optimize the sequences of the candidate peptides and reduce the 

1 0 probability of generating nonfunctional peptide ligands. 

In addition to the use of the twenty naturally occurring amino acids, other 
potential peptide or peptide analogue molecule elements may be used that can be put in 
relationship to the sequence patterns of physical properties determined using the methods 
indicated. For example, D-amino acids or modified and/or pseudo amino acids (e.g., 

15 amino acids bearing acetyl, glycosyl, thiol, chorine, flourine, bromine, alkoxyl, amino 
alkyl, or sulfoximine groups, further including those that are alkylated, acylated, 
methylated and further including those pseudoamino acids that are polycarbonate, 
polyesters, phosphinic, cyclic and others with peptide bonds replaced by a variety of 
other linkages) may be included in the pool of components used to generate the candidate 

20 peptides. Furthermore, other, non-naturally occurring amino acids, dipeptides, 

tripeptides, and the like may be used, as well as non-amino acid compounds. Examples 
of the non-naturally occurring amino acids include, for example, anserine, citrulline, 
cystathionine, homocysteine, 5-hydroxylysine, hydroxyproline, methylhistidine, 
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norleucine, ornithine, phosphoserine, sarcosine, taurine, hypotaurine and other rare amino 
acids. In addition, compounds that involve non-peptide bonds between the constituents 
may be employed if they produce a desirable result, such as increased stability, resistance 
to proteolysis, or increased binding, modulation, activation and/or inhibition of the target 
polypeptide. The only requirements for use of these amino acids and non-amino acid 
components in a manner similar to that of the twenty naturally occurring amino acids in 
the methods of the present invention are that, first, incorporation of the modified amino 
acids and/or components into a linear amino acid chain must be possible, and second, that 
the values for the free energy of transfer of the components (or other of the above listed 
and possible ordered physical properties) must be computable, have quantitatively 
orderable properties relative to one another and be consonant with their assignment as 
dictated by the sequential pattern descriptors such as eigenvector weighting partitions 
such that the component may be assigned to its proper physicochemical group. 

The present invention is illustrated in terms of the following examples, which are 
intended to be descriptive only and is not intended to limit the invention in any way. 



Example 1: 

The 443-amino acid long isoform of the human dopamine D 2 (D 2 DA) receptor 
was transformed into a real numbered AGh P series, H h using the Eyring-Tanford 
hydrophobic^ scale. This H t series (and its all poles maximum entropy power spectral 
transformation, S(co), see below) demonstrated a multimodal distribution (Fig. 2A). In 
place of the a priori selection of orthonormal transformations such as Fourier or Bessel 
functions with which to decompose the receptor's H u i = 1 ■ -443, orthogonal functions 
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were generated from the receptor's Hi directly using the Broomhead-King ("B-K") 
decomposition derivative of methods often named after Karhuenen and Loeve ("K-L"). 
A K-L decomposition of the Hi series of the D2DA receptor involves the autocorrelation 
matrix, Ag 9 of the entire Hu i - 1 • . .443 series, yielding an eigenvector template for D 2 DA 
5 targeted peptides as long as the receptor itself. In the B-K procedure, the H t sequences 
were used to generate an empirically chosen M-lagged data matrix, from which MxM 
covariance matrices, Cm, were computed and decomposed into sets of / orthogonal 
eigenfunctions, where / = 1 . . MJ = 1 . . M. As seen below, this linear 
decomposition yielded eigenvector templates for amino acid assignment of length M. 

10 From the lagged data vectors, and where k = iV-Af+1, the sequence-averaged 

dyadic product, { Hi Hi } > was used to obtain the autoco variance matrix, a M x M matrix, 
C M = l/k{ HiH T i }, using M= 15. We computed the ordered eigenvalues, {v/} and 
the associated eigenvectors, Xfy), of Cm, where i = 1 . . . M and labels the eigenvector, and 
j = 1 . . M refers to the jth component of the eigenvector Xu The eigenvalues, { v J ^um> 

1 5 were ordered from largest to smallest and constituted the eigenvalue spectrum of Cm* 
The similarly ordered and associated eigenvectors, XjQ% were convolved with Hu 
H 2 , .-,H N generating where /= 1...M labels the eigenvector and the j = 1 . JV-Af+1 
(or j=l . . .Abusing the alternate computational form of T^j)) indexed the eigenfunction's 
jth component. The convolution of each of the leading eigenvectors with the Hi series 

20 was performed by computing the sums of the scalar products of the M-length eigenvector 
with an M-length of the Hi series to produce a point in the eigenfunction. Similarly, we 
can sum the scalar products of the eigenvector and a point in the Hi series, giving our 
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alternate computation. Either process was translated down the H t series by one step and 
repeated to generate each of the sequential points of the eigenfunction that corresponds to 
its ordered eigenvalue-associated eigenvector in the computation. We have found that 
when Af » 15, the least squares error was minimized in a fit of the leading eigenfunction, 
5 Yi, dominated by the D2DA receptor's hydrophobic TMs, to the n-block averaged pattern 
of hydrophobic variation, usually called the hydropathy plot. This leading eigenfunction 
demonstrated approximately seven transmembrane segments, and its all poles maximum 
entropy power spectral transformation (iS(co)) demonstrated an average amino acid 
wavelength peak of > 50 amino acids (Fig. 2B). A data matrix of Af * 15 also contained 

1 0 sufficient information such that the secondary D 2 D A receptor eigenfunction, ¥2, could be 
determined to exhibit two putative receptor AGh P binding/modulating mode or modes of 
8.12 and 2.61 amino acids, as seen in its S{ri) (Fig. 2C). The eigenvector associated with 
the secondary eigenfunction, Xz 9 demonstrated all poles, maximum entropy power 
spectra, £(©), with putative D2DA receptor AGhp bindingfaiodulating modes of 8.16 and 

15 2.67 amino acids, as seen in its S(<o) (Fig. ID). These binding/modulating modes were 
closely matched with the modes of the D 2 DA receptor native peptide ligands, such as 
neurotensin, which has an S(o) peak of » 8.13 amino acids. M= 15 is within the middle 
of the « 5-30 amino acid length range of most physiologically active peptides. Most 
peptides with the capacity to bind antibodies and elicit an antibody response are also in 

20 the range of about 5-30 amino acids in length. 

Figures 3 A and 3B are two-dimensional graphical representations of the Morlet 
wavelet W(a,b) transformation of the Hi of the D 2 DA receptor. In these graphs, sequence 
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position is graphed along the x-axis, phase amplitudes along the y-axis and m =fldd) is 
fixed at the two characteristic peaks (hydrophobic free energy binding/modulating mode 
or modes) of the S(q) transformation of Y 2 , as well as at the highest phase amplitudes of 
the W{a y b) transformations of the Hi of the D 2 DA receptor, at tn, co » 2.3 and 8.1 amino 
5 acid residues. Figures 2 A and 2B demonstrate that although both the 2.3 amino acid and 
the 8.1 amino acid wavelengths of the D 2 DA receptor have phase amplitude peaks that 
are distributed throughout the Hi length of D 2 DA, the most prominent of the 8.1 amino 
acid phase amplitude sequence locations (marked by arrows) correspond to the 
extracellular loops EL-I, between TM 2 and TM 3 (« residues 85-105); EL-II, between TM4 

10 and TM 5 ((« residues 190-210); and EL-III, between TM 6 and TM 7 (« residues 390-410). 
The brain peptide neurotensin is believed to mediate its actions through the D 2 DA 
receptor, and neurotensin exhibits an S((o) peak of co" 1 » 8.13 amino acids, which matches 
well with that of the D 2 DA receptor. 

Peptide construction from the eigenvector template derived from the D 2 DA 

1 5 receptor was performed with the y-axis of X 2 as graphed in Figure 2D (left) being divided 
into four equal segments corresponding to the natural 4-partition of the AGh P values of the 
twenty naturally occurring essential amino acids listed above in Tables 1 and 2. 
Probability weightings for amino acid members of each of the four AGh P groups were 
assigned on the basis of their relative occurrences in human cerebrospinal fluid (CSF), 

20 reflecting the brain's amino acid pool available for peptide synthesis. In addition, 

probability weightings were assigned on the basis of the amino acid distribution in each 
of the four groups of neurotensin, which we have shown previously to modulate the 
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kinetics of binding by the human D 2 DA receptor. Based on these distributions, weighted 
random assignment of amino acids to each of the 15 points of the 4-equipartitioned X 2 
generated the new peptides. The first two peptides were derived from the CSF pool 
probabilities, SHQRWEYKGVNCPVY ("SHQR"; SEQ ID NO:l) and 
THQ AFHYCNKQCLVI ("THQA"; SEQ ID NO:2) (Table 3), and were synthesized to > 
95% purity (as determined by HPLC and mass spectrometry) by Multiple Peptide 
Systems (La Jolla, C A). Two additional peptides using an idealized X 2 and with 
probability weightings derived from the amino acid composition of neurotensin rather 
than human CSF, ERNRKPLRPKNKYLI ("E. . .PL"; SEQ ID NO:3) and 
ERNRKPYRPKNKYLL ("E. . .PY"; SEQ ID NO:4) (Table 3), were also designed and 
synthesized for microphysiometric testing. The last eight D 2 DA targeted algorithmically- 
derived peptides were produced using the X 2 eigenvector of the M = 1 5 covariance 
matrix, C M , of the human, long isoform, D 2 D A receptor as the template for amino acid 
assignment. 

As an example of one of many possible physiological assays that may be used to 
evaluate the actions and potencies of designed peptides, two independently derived cell 
systems were examined with respect to the peptide action and/or modulation of their 
external acidification rate ("EAR") to dopamine. The mouse LtK fibroblastoma cell 
system was generously provided by Frederick Monsma (Hoffman-LaRoche, Basil, 
Switzerland). The CHO (Chinese hamster ovary) cell system was generously provided 
by Richard Mailman (Univ. of North Carolina, Chapel Hill, NC). Both cell systems were 
stably transfected with human long isoform D 2 DA receptor cDNA, which had been 
isolated from a human striatal cDNA library, sequenced and subcloned into the 
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expression vector pRC/RSV (Invitrogen). The transformed Ltk system was characterized 
by lower baseline responsivity to its native agonist, dopamine, as measured in total milli- 
pH units (mpH), In contrast, the transformed CHO system manifested a higher baseline 
responsiveness to dopamine. Both systems were grown to confluence in DMEM 
5 containing 10 % FBS. The cells were serum-starved 18-24 hours prior to use, and then 
assayed for EAR using a microphysiometer (Cytosensor; Molecular Devices, Sunnyvale, 
CA) in low buffering DMEM with 0.1% culture grade BSA. 

The determination of EAR by microphysiometry involves a proton-sensitive 
silicon semiconductor photocurrent-driven sensor which measures changes in EAR 

10 resulting from effector-evoked alterations in cellular glycolytic and respiratory energy 
metabolism and/or alterations in sodium-hydrogen exchanges across cellular membranes. 
Protonic H*, generated by such energy metabolism or exchanges, neutralizes the charge 
on the surface of the semiconductor, reducing the photocurrent produced at a rate linearly 
related to H + production. 

1 5 The microphysiometer monitors pH in flow-through chambers containing the 

receptor-transfected cells. Generally, if the cells lines used are adherent cell lines, the 
cells are seeded into "capsule cups". If the cell lines are non-adherent cell lines, then the 
cells are immobilized in a fibrin matrix. For all microphysiometer runs, modified low 
buffering DMEM containing 0.1 % BSA is pumped across the cells at a rate of 

20 approximately 100 |il/min s during which time the pH of the microenvironment 
surrounding the sensor surface is maintained at a relatively constant value. The 
measurement of the acid output rate of the cells, termed the acidification rate, is made 
when the fluid flow is periodically halted to allow buildup of acidic metabolites in the 
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chamber, resulting in an alteration in the pH of the fluid. The pH is measured in 
millivolts, and converted to milli-pH units. The changes in pH are expressed as changes 
in milli-pH units per minute following the linear, time-dependent buildup of H + during 
intermittent periods of pump arrest followed by washout. Integration of the EARs over 
5 the time of action of dopamine yields an estimate of the total milli-pH units (measured as 
the area under the curve by trapezoidal approximation) generated during the action of the 
natural ligand alone, compared with that of the ligand when preceded by the infusion of 
the algorithmic peptide. This data is plotted as average sensitivity in the range of 0.001 
pH units, and changes as little as 2 % of the control are reproducibly detectable. Ligand 

10 induced, receptor-mediated increases in cell metabolic and Na + -H + membrane regulatory 
activity is seen as an increase in the acidification rate. 

Dopamine was infused at concentrations approximating its EC50 in this system, 
that is, « 1 |aM. Following pilot studies which indicated a consistency in sensitivity and 
direction of effect, the twelve peptides were surveyed at 1 \xM concentrations. Small 

1 5 Kolmogorov-Smirnov distances supported the assumption of normality in all of the data 
sets, so within chamber-paired, one-tailed t-tests with a significance criterion of p = 0.05 
were used. 

Figures 4A-4D summarize the EAR responses to dopamine infusion with respect 
to the influence of SHQR and THQA in the two D 2 DA receptor-transfected cell systems, 
20 in which the former significantly potentiated the dopamine-induced increment in total 

milli-pH units in both cell systems. We report the results of one-tailed t-tests with pairing 
within chamber as t(#>, where # represents the degrees of freedom of the paired 
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comparison and p denotes the probability of such results occurring by chance. For the 
SHQR peptide in the LtK system, t (3) = 13.28, p - 0.0009, and for the SHQR peptide in 
the CHO cell system, t (3) = 28.06, p < 0.0001. THQA did not significantly potentiate the 
dopamine response in either system, t@) = 0.620 and tp) = 1.309, p > 0.05, respectively. 
5 Figures 5A-5D contain graphs of the influence of the peptides E. . .PL and E, . .PY on the 
EAR response to dopamine in the two D 2 DA receptor-transfected cell systems. Both 
peptides demonstrated statistically significant activation, t(7)~ 25.47, p <0.0001 and t@) = 
69.830, p < 0.0001, respectively, in the LtK system. However, neither of the E. . .PL and 
E...PY peptides influenced the dopamine-induced EAR of the CHO cells significantly, 

10 with t(3) = 1.542, p > 0.05 and t (7) = 1.283, p > 0.05, respectively. Three of the 

remaining eight peptides exhibited statistically significant effects on at least one of the 
two receptor-transfected cell systems (Table 3). The overall "hit rate", as measured by 
modulation of the kinetics of the EAR of two transfected cell lines to dopamine, for these 
peptides was thus 50% (i.e., six of twelve peptide candidates that were synthesized and 

1 5 tested statistically significantly altered EAR in one or both of the D 2 DA receptor- 
transfected cell systems used). All D 2 DA targeted peptides whose effects reached 
significance increased EAR. 

A set of EAR dose response curves were computed for SEQ ID NO:l across 
concentrations of dopamine (10' 8 5 M to 10" 55 M) and the peptide SEQ ID NO:l (10 nM 

20 to 3 jaM) (not shown). LtK cells were used for these experiments. The resulting dose 
response curves manifested asymptotic sigmoidal kinetics, suggestive of positive 
cooperativity. 
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Tables 3, 4 and 5 show the sequences of the various peptides synthesized by the 
methods of the present invention and their effect on the cell test systems. 



TABLE 3 



HT IM AN DOPAMINE D? (D?DA) RECEPTOR TARGETS 


D PEPTIDES 


SEQUENCE 


SEQ ID 
NO 


DIR. 
EFF1 


2CT 


MODULATORY 
EFFECT 


CHO 


LtK 


CHO 


LtK 


H-SHQRWEYKGVNCIVY-OH 


1 






**# 


*** 


H-THQAFHYCNKQCLVI-OH 


2 


o 

I 




ns 


ns 


H-ERNRKPYRPKNKYLL-OH 


3 


? 


o 

I 


ns 


s|s s|e s#e 


H-ERN KLIN Y IsJN isJN Js. Y lL-KJn. 


A 
*t 


? 


? 


ns 


*** 


H-SHTAYHWMSCGKTVI-OH 


5 


ns 


ns 




* 


H-SRQAFHYKNVQVLVL-OH 


6 


? 


? 


ns 


ns 


H-SHQAWRYKNVNCYVI-OH 


7 


? 


ns 


? 


*** 


H-GETAFRYVNCNVYVYI-OH 


8 


** 


*** 


ns 


ns 


H-GHSAWRWKSKNVYMI-OH 


9 


ns 


ns 


ns 


ns 


H-NASALHLVGVQCWVY-OH 


10 


? 


ns 


? 


ns 


H-SWQAIRICQKGVLMY-OH 


11 


? 


ns 


? 


ns 


H-SHSRWRTVSVNVLCY-OH 


12 


? 


ns 


9 


ns 



*0.01<p^0.05, **0.001<p<;0.01, ***p<0.001. ?=not yet tested, ns=not significant. 



Example 2: 

Peptides derived from receptor protein systems other than the D 2 DA receptor 
1 0 were also tested for their effects on their respective receptors. For the human muscarinic 
Ml receptor, CHO cells were transfected with the muscarinic Ml receptor cDNA derived 
from a human cDNA library essentially as described by Buckley et al. (Mol. Pharmacol. 
1989 35:469-476). Briefly, the coding region of the Ml receptor was obtained from a 
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human cDNA library and cloned into the expression vector pcDNA3 (Invitrogen, San 
Diego, CA). CHO-K1 cells were transformed with the construct, using the calcium 
phosphate method. Stably expressing transformants were obtained in the presence of 250 
jig/ml geneticin. Transformed cell lines expressing the human NGF receptor also were 

5 obtained. The effects of the peptides derived by the methods of the present invention on 
the activities of the corresponding receptors in the transformed cell lines were evaluated 
in the same manner as described above for the D 2 D A -targeted algorithmically derived 
peptides, using the EAR test system. 

In the case of the Ml receptor, ten peptides were obtained, using the methods of 

1 0 the present invention. Of these ten peptides, five (50%) had a statistically significant 
effect on the EAR due to carbachol in the Ml receptor-transfected CHOK1 cells (e.g., 
Figs. 6A-6B, Table 4). All these effects were direct or modulatory decreases in EAR. 
This contrasts with the positive direct or modulatory effects of the tested peptides on the 
EAR to dopamine in D 2 DA receptor-transfected cell lines. 
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TABLE 4 



HITMAN MUSCARINIC Ml RECEPTOR TARGETED PEPTIDES 


SEQUENCE 


SEQID 
NO 


DIRECT 
EFFECT 


MODULATORY 
EFFECT 


H-FSFQCKSESfYEALGY-OH 


1 "? 
13 


** 


** 


H-FSFGVKSWQYHALGY-OH 


1 A 

14 


ns 




H-ITFTVKGLTLAAFTY-OH 


15 


o 

f 


*** 


H-ISFNKCTWSFERYSL-OH 


16 


ns 


* 


u T7MT CVVnWMYT? A YNT -OH 


17 


ns 


** 


H-LNYQKKQYTYAAWQF-OH 


18 


ns 


ns 


H-LTYGVMNYGFAAFGF-OH 


19 


ns 


ns 


H-LGFSVCPITLAELTY-OH 


20 


ns 


ns 


H-LGLGVCPINLAALTW-OH 


21 


? 


ns 


H-LTWNVKTYSLHELPL-OH 


22 


ns 


ns 



*0.01<p<0.05, **0.001<p<0.01, ***p^0.001. ?=not yet tested, ns=not significant. 



For the NGF receptor, 1 1 peptides were obtained, using the methods of the 
present invention. Of these 1 1 peptides, eight (73%) exhibited a statistically significant 
change in EAR (Table 5). 
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TABLE 5 



HUMAN NERVE GROWTH FACT 


fOR RECEPTOR TARGETED PEPTIDES 


SEQUENCE 


NO 


EFFECT 


MODULATORY 
EFFECT 


H-DLCRSARSDIEVTEY-OH 


23 


** 


*** 


H- RFVASAATEIEVNRL -OH 


24 


ns 


♦ * 


H-HYCASADPRIHKNAL-OH 


25 


ns 




H-DFVDGAAGRLHKGEY-OH 


26 


ns 


** 


H-DEKATEATDIEKGHL-OH 


27 


ns 




H-RFVDNDATDIEKGRI-OH 


28 


* 




H-RFVRGDRNHFDCGEL-OH 


29 


* 


*** 


H- HFVRNERTHFDVSAL -OH 


30 






H-AYKHNEATDIEKGDF-OH 


31 


ns 


ns 


H-HKRKEATHIEKSAL-OH 


32 


ns 


ns 


H- HIVEGRAPEIACGEY -OH 


33 


ns 


ns 



*0.01<p<0.05, **0.001<p<0.01, ***p< 0.001. ?=not yet tested, ns=not significant. 



Thus, 33 total peptides obtained using the methods of the present invention, for all 
of the receptor systems tested. Of these, 19 had a significant effect on the EAR of the 
transformed cell lines directly or in response to the native ligand, resulting in an overall 
hit rate of 57.6%. At a rate of 5 per 100,000,^(5) = 0.00005, as the random 
combinatorial prior probability of hits, and 2 per 4, p(A) = 0.5 as the probability of 
physiological action observed of eigenvector template-generated peptides, a Bayesian 
theorem says that the latter would occur under conditions of the former like: 

pjA 1 B)p(B) = 0.000025 x 0.00005 _q 25-1Q" 8 . 
p{A) 0.5 
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Thus, an overall average hit rate of 57.6% achieved by the receptor-targeted 
algorithmically-derived peptides produced by the methods of the present invention 
appears to be orders of magnitude more efficient for lead peptide generation when 
compared to the conventional methods of randomly generated peptide libraries. 

5 Cross-over experiments were performed to determine the specificity of the active 

peptides for the receptor protein from which they were derived. When the D 2 DA targeted 
algorithmically-derived peptides that had a significant effect on EAR to dopamine in 
D 2 DA receptor-transfected cell lines in were tested for their influence on the EAR to 
carbachol in Ml receptor-transfected CHO-K1 cells, no effect was observed. Similarly, 

1 0 no effect on the EAR to dopamine in D 2 D A receptor-transfected cell lines was observed 
in the presence of the peptides that exhibited a negative allosteric effect in the Ml 
receptor-transfected cell lines. Therefore, the peptides appear to be selective for the 
mode-matching receptor proteins from which they are derived. 



15 Example 3: 

Using the redundant subsequence template method described above, peptides 
were derived from the known polypeptide calcitonin. The parent family of known 
calcitonins are 31 amino acids in length, which was reduced to 10 amino acids using the 
redundant subsequence template method to produce the peptides listed in Table 7. The 

20 redundant subsequences were generated by examination of the calcitonin sequences of 
eight different species (Table 6). 
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Table 6 





Fnnr-nnmher Hvdronhobic Free Enersv 
Codes 


Nonoverlapping Repeated 
Subsequences 


Human 


3113113331141124134214411241312 


4112413; 311 


Swine 


3113113331244213114224113141421 


1131; 421 


Cow 


3113113331244323114224113141421 


1131; 142 


Sheep 


3113113331244323114224113141421 


1131; 142 


Rat 


3113113331141123134214411141312 


1141; 3112 


Eel 


3113113331331123233114421231211 


3311; 311 


Salmon 


3113113331331123233114421111111 


3311; 311; 111 



Conventionally, calcitonin is administered by daily injections to post-menopausal 
women suffering from osteoporosis. By reducing the peptide length from 31 to 10 amino 
acids, the resulting peptides are more easily administered by transdermal and inhalation 
methods. The peptides listed in Table 7 are examples of peptides generated from a 
human non-overlapping redundant subsequence template (i.e., 31 141 12413), and are 
weighted by the amino acid distribution of the human calcitonin receptor. 
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Table 7 



Examples of Human Calcitonin-Targeted Peptides from Non-Overlapping Redundant 
Substring Template of Human Calcitonin 


VPXTT "DXTET XTfcT 
JSxiN Lr IN lS.LiN rv 


V 1 IN W IN VJIvU > J\. 


VOTYPPHFPV 

Vy l ijrxj.jj.jr v 


KTTINGHISK 


V xiNlAjiNxjltjrV 


MOWFPT A TMV 


KGNT NTDT NM 


VGGYGTDYNM 


CiMMrorDll V 


Vr alv^vjrl i vjM 


VTPT ^ATXTK' 


MOGYTNDTPV 

lYlV^VJ I 111 J-/ JX V 


JVL 1 1 1 rir 


VVjxs.li vri i i Jv 


V IN IN lvJjC I IN V 


VNOWONHYTM 

V 1>V^ TV Yilll X X XTX 


VN lrOlbLoC 


Vrrf llNrlVVV^Jv 


MPPWP^nYPf 
ivir i W i i r 


KPTFSNAYNV 


CNNIGNRLSC 


VGTLNPAFSV 


KOSFQSELNK 


VTNFSNALSM 


KGNFTPEWPC 


CGNYGTRFSK 


VPSLTTRLQV 


VTPINSEFPC 


MGPLPQAFQC 


CSSLQOALTV 


VQPLQGHLPV 


KNOLNTHIGK 


KSNIGPALTM 


MPSIPTHLNK 


VSOFNQAWGV 


VOSINNAIGK 


VSQYGQELQV 


KNNYGOAFTV 


VPSLNSALGV 


MGTFQPDWQV 


VSPYQSHFNV 


KNQLNTEINC 


MNSIQTDFTM 


VQTISSRWGK 


MGGWGPALNC 


KNPLNNHLNM 


VOSLTNDISK 


MGNITQDLQC 


CTGYTNAIQM 


VNGIGQAINV 


KGNINPAYNV 


KGSYTTELGV 


MNTLQQAYPK 


CPGITGDFQK 


KTGLNNEINV 


KNSYSPELTV 


VOPYNGELNM 


MTOFQSHITV 


VOSFTNEIQC 


CNSYTPEFPC 



The hydrophobic wavelength of a peptide containing L- amino acids is the same 
5 as a peptide containing D-amino acids in which the sequence is inverted. Such "retro- 
inverso" peptides have been previously described (Chorev, M. and Goodman, M. (1995) 
Trends Biotechnol. 13:438-445), but their use as mode-matched binding peptides has 
never before been contemplated. A retro-inverso peptide containing D-amino acids and 
having the sequence LHGKEIDTAETAKID was synthesized (SEQ ID NO: 28). This 
1 0 sequence of this peptide is inverted from that of peptide of SEQ ID NO:27, used in the 
NGF receptor inhibition assay. The peptide of SEQ ID NO:27 significantly down- 
regulated the EAR response of transformed cell line containing the NGF receptor, at a 
significance level of p < 0.001 . The peptide of SEQ ID NO: 28 was tested in the same 
assay as the L-amino acid, forward sequence peptide of SEQ ID NO: 27. As shown in 
15 Fig. 7, this peptide also down-regulated the EAR response of PC- 1 2 cells to NGF, to an 
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extent comparable to that seen for the L-amino acid, forward sequence peptide of SEQ ID 
NO: 27. . 

Retro-inverso versions of peptide antigens are known to evoke more powerful 
antibody responses than L-amino acid forward sequence versions and the antibody 
5 responses also lasts longer. This provides additional support for the idea that it is the 
hydrophobic mode patterns that largely dictate binding, because the orientation of the 
retro-inverso peptide backbones are completely altered with respect to that of the forward 
sequence peptides, but the hydrophobic mode patterns are not. As a result, "hydrophobic 
mode matched" retro-inverso peptide antigens could be designed to have stronger 

10 immunogenic properties than the usual peptide fragment of proteins used as antigens. 
Such retro-inverso peptide antigens could be orally administered, since they would be 
resist proteolyetic digestion. However, there is also the possibility of a patient 
developing resistance to their effects due to the generation of antibodies against such 
peptides. Such a response may differ from one retro-inverso, hydrophobic free energy 

1 5 mode matched peptide to another. 

The methods of the present invention may be used to produce peptides useful in 
variety of investigative, therapeutic and diagnostic applications as listed in part in the 
examples listed above. In addition to these applications, the peptides may be used in the 
detection and/or treatment of cancerous tumors. The peptides may also be used in the 

20 detection and/or treatment of various other disease conditions, and may also be useful in 
the detection of contaminants in food, water or soil It will be appreciated that if the 
sequence of a particular polypeptide that is specifically or exclusively associated with the 
disease condition, tumor, or contaminant is known, then peptides that will bind, modulate 
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the function of, activate or inhibit those polypeptides may be synthesized by using the 
methods of the present invention. When used to treat a tumor, the peptide may be 
conjugated to or incorporate a cytotoxic agent, such as a radioisotope or a toxin. When 
used for detection, the peptides may be conjugated to a molecule that can be visualized or 
5 otherwise detected, such as a radioisotope, a chromophore or a fluorophore. The peptides 
of the present invention may be used to screen bodily samples for the presence or absence 
of a particular polypeptide. Examples of such bodily samples include blood, plasma, 
blood products, urine samples, fecal samples, tissue biopsy samples, skin samples, semen 
samples, and epithelial cell samples. When used to screen for tumors or disease 

10 conditions, or when used as a therapeutic, the peptides may be included as a component 
in a diagnostic or therapeutic kit, respectively. The peptides may also be used in areas of 
research, such as molecular biology, pharmacology, neurobiology, intracellular signaling 
and the like, to explore the functions and pharmacological responsivities of proteins, 
polypeptides or peptides of unknown functions. For example, a tissue culture cell line 

1 5 transfected with a cloned orphan receptor may be incubated with various mode-matched 
peptides and tested for any number of cellular activities that may be associated with that 
receptor. The use of peptides in general in such applications are well known to those in 
the art; therefore the peptides produced by the methods of the present invention may be 
used in the above-cited applications in the usual manner, without the need for undue 

20 experimentation. 

Although the invention herein has been described with reference to particular 
embodiments, it is to be understood that these embodiments are merely illustrative of 
various aspects of the three template generating methods and the amino acid assignment 
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methods of the invention. Thus, it is to be understood that numerous modifications may 
be made in the illustrative embodiments and other arrangements may be devised without 
departing from the spirit and scope of the invention. Throughout this application various 
publications may be cited. Where cited, the contents of these publications are hereby 
incorporated by reference into the present application. 
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